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Abstract 

We present the Callan-Symanzik-Lifshitz method to approaching the critical behaviors of systems 
with arbitrary competing interactions. Every distinct competition subspace in the anisotropic cases 
define an independent set of renormalized vertex parts via normalization conditions with nonvan- 
ishing distinct masses at zero external momenta. Otherwise, only one mass scale is required in the 
isotropic behaviors. At the critical dimension, we prove: i) the existence of the Callan-Symanzik- 
Lifshitz equations and ii ) the multiplicative renormalizability of the vertex functions using the in- 
ductive method. Away from the critical dimension, we utilize the orthogonal approximation to 
compute higher loop Feynman integrals, anisotropic as well as isotropic, necessary to get the expo- 
nents rjn and Vn o,t least up to two-loop level. Moreover, we calculate the latter exactly for isotropic 
behaviors at the same perturbative order. Similarly to the computation in the massless formalism, 
the orthogonal approximation is found to be exact at one-loop order. The outcome for all criti- 
cal exponents matches exactly with those computed using the zero mass field-theoretic description 
renormalized at nonvanishing external momenta. 
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I. INTRODUCTION 



Universality classes characterizing the ordinary critical behavior of systems undergoing 
phase transitions jl| have an interesting parallel with those associated to competing systems 
of the Lifshitz type The ordinary universality hypothesis states that all universal 

amounts characterizing the transition like critical exponents, amplitude ratios of certain 
thermodynamic potentials above and below the critical temperature, etc., do not depend on 
the microscopic details of the systems but depend upon the number of components of the 
order parameter and the space dimension d of the system. Lifshitz universality classes, on the 
other hand, depend on an additional parameter : the number m of spatial directions where 
competition takes place. The simplest competing systems belonging to the m-axial Lifshitz 
universality classes correspond to m = 1 (uniaxial) and are examples of complex systems 
with two ordered phases as well as one disordered phase in the vicinity of the Lifshitz critical 
point. The m space directions are called the competition axes. Turning off the competing 
interactions is equivalent to taking the limit m — and the competing system should turn 
into an ordinary critical system. 

Lifshitz multicritical points appear at the confluence of a disordered phase, a uniformly 
ordered phase and a modulated ordered phase. The spatially modulated phase is character- 
ized by a fixed equilibrium wavevector, which goes continuously to zero as the Lifshitz point 
is approached. But how can we realize these ideas in terms of concrete models describing 
actual physical systems? The language of magnetic systems is particularly convenient to 
find a simple realization of this critical behavior in terms of a lattice model named ANNNI 

n 

model It is a d-dimensional Ising model with ferromagnetic couplings between first 
neighbors together with antiferromagnetic exchange forces between second neighbors along 
only one space direction. Competition in the interactions arises by varying their strengths. 
More precisely, the Lifshitz point corresponds to a particular value of the ratio between the 
ferromagnetic and antiferromagnetic exchange interactions. This situation corresponds to 
the anisotropic uniaxial criticality (m = 1). The ANNNI model can be applied to describe 
the critical behavior of many systems. However, experiments js, 6, 7, ^ as well as theoretical 
investigations 9|] have confirmed that the magnetic material MnP yields a simple realiza- 
tion of the three-dimensional ANNNI model: it has a pure uniaxial (A^ = = 3,m = 1) 
Lifshitz critical behavior. 
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The competing axis can be classified according to the number of neighbors which are 
coupled via competing interactions. The ANNNI model can be generalized to include m 
space directions with competing interactions. We employ the notation 1712 = m to label this 
subspace according to the number of neighbors they connect through competing exchange 
forces. In that case the wavevector characterizing the modulated phase has m2 components. 
The critical system under study presents an m2-fold (or m2-axial) Lifshitz critical behavior, 
which can be either anisotropic m2 < d or isotropic d = 1712 (close to 8). 

If the field (order parameter) has components, the Lifshitz universality classes are 
defined by the triplet {N,d,m2). These universalities correctly reduce to the Ising-like 



universality class {N,d) in the limit m2 lOj. These criticalities have encountered 



applications in many real physical systems like liquid crystals 



polymers 13|, microemulsions 



[14 ], high- 



lll . ferroelectrics [12] , especial 



Tc superconductors [15!], magnetic materials 



71, l8|, etc. Furthermore, other aspects have been studied like the formulation of quantum 
phase transitions in Lifshitz points as well as the connection of Lifshitz type field 

theories with weighted scale invariant quantum field theories 19]. 

From the technical point of view, the original calculation of usual critical exponents 
belonging to the Ising-like universality class (without competition) were performed using 
the renormalization group and e-expansion methods via diagrammatic perturbation of field- 
theoretic renormalized massive theories ij]. Some time later this method was reformulated 
such that the former approach was reduced to the computation of a few diagrams (IPI ver- 
tex parts) yielding the same exponents through the use of a renormalized massless scalar 
field theory at nonvanishing external momenta 



21 



221 ]. Inspired in this massless frame- 



work, many ca' 



forward 



24 



culational schemes in the study of m-axial Lifshitz points have been put 



251 ]. In particular, the unconventional renormalization group arguments 



in the anisotropic criticalities along with analytical solution methods to resolve Feynman 
path integrals have permitted a better comprehension of the critical properties of Lifshitz 
points using the massless scalar field-theoretic setting renormalized at nonvanishing external 
momenta together with the renormalization group equations in the large distance infrared 



regime 



26, 



27]. 



The method of Refs 



21 



221 ] to obtain critical exponents of ordinary critical systems 
was adapted to include the massive theory renormalized at zero external momenta. It al- 
lowed to understand the connection between the infrared behavior of the solution to the 
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renormalization group equation in the massless theory with the ultraviolet behavior of the 



solution to the Callan-Symanzik equation 
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2^ in the massive theory 30|] . These calcula- 



tions are more involved, since the massive Feynman integrals are intrinsically more difficult 
to solve than their counterparts in the massless framework. The benefits of this study is that 
a more direct connection with quantum field theory can be made in the ultraviolet regime. 
Besides, this extra information is useful to have a proper understanding of other universal 
quantities like the equation of state, amplitude ratios of certain thermodynamic potentials 
above and below the critical temperature, etc. We have shown in a previous work 3l|] that 
it is possible to formulate an appropriate renormalized field-theoretic setting for massive 
scalar fields with quartic self interactions in order to compute the critical exponents asso- 
ciated to m-fold Lifshitz points. There were two main steps explicit in that construction. 
The ffist one is the appearance of two independent mass parameters necessary to describe 
the two inequivalent space (or momentum) directions present in anisotropic criticalities in 
an independent manner. Consequently, the Callan-Symanzik-Lifshitz equations allowed to 
solve the problem at the repulsive ultraviolet fixed points in the anisotropic cases with two 
independent mass scales. Second, the isotropic universality class only needs one mass scale. 
Needless to say, the results for the exponents using either the massless approach or the 



massive method are the same 



3l|. 



A different generalization of the ANNNI model can be considered to include further alter- 
nate competing interactions, for instance, up to third neighbors. From a phenomenological 
viewpoint, the ffist nontrivial example of a higher character Lifshitz point (see below) occurs 
for a uniaxial third character Lifshitz point. The three-dimensional phase diagram consists 
of two parameters varying with the temperature, i.e., the ratio of exchange interactions be- 
tween the second and nearest neighbors as well as the ratio of exchange couplings between 
the third and the second neighbors. When the temperature axis is projected on the plane 
of these parameters, there is a region of intersection where the different phases associated 
to the system encounter each other in the uniaxial Lifshitz point of third character, whose 
existence was established numerically 32| • As before, we should emphasize that the uniaxial 
third character competing axes connect up to third neighbors with alternate exchange forces. 
When these space dimensions occur along directions, the critical competing system is 
said to represent the ms-axial third character Lifshitz critical behavior. These universality 
classes are characterized by (A^, rf, ma). The uniaxial third character Lifshitz criticality is 
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the particular case of the ms-axial behavior for ma = 1. 

Using these ideas, the generahzation of competing systems whose competing axes hnk 
up to L neighbors via alternate exchange forces can be easily understood. Consider a d- 
dimensional Ising model with competition interactions connecting L neighbors. Let mi 
be the competition spatial subspace where exchange short range couplings take place with 
ferromagnetic interactions between first neighbors, antiferromagnetic forces between second 
neighbors, ferromagnetic interactions between third neighbors and so on, with alternating 
signs for the exchange forces up to the Lth neighbors. This universality class is characterized 
by {N,d,mL) [33|, |34]- Even though the anisotropic {rriL < d) and isotropic behaviors {d = 
iriL) can still be defined, the anisotropic situation can be described with only 2 independent 
scales. Thus, this model does not correspond to the most general competing system. 

Generic competing systems of the Lifshitz type have been introduced recently. A simple 
realization of their critical properties using the lan guage o f magnetic systems can be visu- 



alized through a lattice model called CECI model [35|, ISGj. It is a generalized Ising model 
with several distinct types of competing axes. Each competition subspace is perpendicu- 
lar to each other. It describes a complex system, which in a simple situation possesses L 
ordered phases, as well as one disordered phase, near the generic higher character Lifshitz 
point. This is so because the CECI model contains simultaneously independent competing 
axes whose exchange couplings are independent in each spatial subspace. Consequently, 
in the anisotropic situation there are m2,ms, ...,mL types of competing axes, as well as 
{d — 7712 — ... — 7711) space directions where there are only ferromagnetic couplings among 
nearest neighbors. We then define mi = d — 7712 — ... — rriL in order to unify the treatment to 
all subspaces, the competing and noncompeting ones. The noncompeting subspace can be 
viewed as the competing axes with only first neighbors interacting ferromagnetically. The 
ordinary critical behavior can be understood in terms of this generic competing system as 
the isotropic particular case for d = mi with m2 = ... = m^ = 0. 

The generic L-th character Lifshitz anisotropic universality classes are now identified by 
the grid (A^, mi, m2, m^) (actually the same as {N,d, 7712, ...,7711)), therefore generaliz- 
ing the L-th character Lifshitz universality class previously discussed, which depends on 

{N, d, 771l) 



33 



34| . However, can we find physical systems which are realizations of higher 



character Lifshitz points, or does this model have only academic interest? For instance, 
higher character Lifshitz points can show up in blends of diblock copolymers. In fact, the 
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earlier nomenclature regarding the higher character Lifshitz points can be translated into 
the modern terminology if we identify the L-th character Lifshitz point with the (former) 
Lifshitz point of order {L — 1) 33|. Thus, it can be verified that Lifshitz points of up to 6th 
character are possible in those blends of diblock copolymers 37|]. This is an encouraging 
evidence that competing systems represented by the CECI model might find applications in 
many other examples of actual physical systems yet to be discovered. The ANNNI model 
can be embedded into the CECI model, such that the former can be retrieved from the 
latter when we switch off the competing interactions beyond third neighbors by taking 
m2 = 1,1713 = ■■■ = ttll = 0. Moreover, the isotropic behaviors have their universality classes 
completely specified by {N,d = rriL), whose critical dimension dc = 4:L depends as well 
on the number of neighbors coupled via competing interactions. They also reduce to the 
isotropic m-axial case when L = 2. Thus, the usual m-axial Lifshitz criticality is a particular 
case of the generic competing system described by the CECI model. 

In this paper we generalize the framework of renormalized massive scalar fields previously 
introduced to compute critical exponents pertaining to the m-axial Lifshitz universality 
classes 3l|] in order to study the most general competing system. It can be mathematically 
understood in terms of an anisotropic field-theoretic description of a massive renormalized 
scalar field theory including arbitrary higher order derivatives. Each higher derivative 
term in the Lagrangian density defines a certain type of competing spatial subspace. There 
are L different types of competing axes which result in L sets of independent masses and 
coupling constants. We derive the Callan-Symanzik-Lifshitz equations with L independent 
mass scales, therefore generalizing the previous massive formulation for anisotropic m-axial 
Lifshitz critical behaviors defined by only two independent mass subspaces jsil. We study 
the solutions of these generalized Callan-Symanzik-Lifshitz equations with several indepen- 
dent mass scales. We show that their ultraviolet behavior at the nonattractive ultraviolet 
fixed point is completely equivalent to the solutions of the renormalization group equations 
formulated in the massless theory with several independent momenta scales at the infrared 
fixed point. We focus on renormalized perturbation theory in order to compute the critical 
indices by diagrammatic means within this technique. The perfect agreement with those 
calculated in the massless framework at the same loop order is a clear evidence that the 
universality hypothesis is obeyed as expected. 

The critical indices rjL and ul associated to each type of competition axes are computed 
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in the anisotropic cases using the orthogonal approximation previously introduced in [351.136 1 
to treat the massless field theory formalism. We utilize several mass scales, one for each com- 
peting subspace, and the corresponding renormalized theory is defined at vanishing external 
momenta. Physically, the masses correspond to the independent correlation lengths ^i,..., 

which go critical simultaneously and naturally take place in the anisotropic criticalities. 

In addition, we use a similar framework to define the isotropic behaviors with only one 
mass scale and derive the corresponding Callan-Symanzik equations for them. We compute 
the critical exponents perturbatively using Feynman graph techniques. The orthogonal 
approximation is employed to compute the analogous critical exponents in the isotropic 
behaviors. Besides, we shall demonstrate that in spite of the highly nontrivial feature of the 
Feynman integrals in the massive theory, they can be computed exactly giving exact results 



35 



36|. 



for the critical exponents identical to those found in 

The normalization conditions are presented in Sec. II. There we motivate the origin of 
the several independent mass scales in the anisotropic behaviors and show how they can be 
understood in terms of the phase diagrams of the CECI model. A simpler treatment will be 
given for the isotropic behaviors as well with only one mass scale for each type of competing 
axes. 

We discuss the one-loop renormalizability at the critical dimension of the various 
anisotropic and isotropic situations in Sec. III. Using these concepts, we prove inductively 
the finiteness of the multiplicatively renormalized vertex parts at all orders in a perturbative 
expansion and demonstrate that the Callan-Symanzik-Lifshitz equations exist in Sec. IV. 

In Sec. V we discuss the Callan-Symanzik-Lifshitz (CSL) equations slightly away from 
the critical dimensions {d = d^ — eij.We derive those equations with several mass scales for 
the anisotropic criticalities. We present the solution of the CSL equations in the ultraviolet 
regime and show that at the ultraviolet nonattractive fixed point it has the same scaling 
form as the solution of the renormalization group equations in the infrared regime. We show 
that the critical exponents calculated by diagrammatic means can be identified with the 
anomalous dimension of the field and that of the composite operator at the ultraviolet fixed 
point. 

The computation of the anisotropic critical exponents are presented in Sec. VI. The results 
for Feynman graphs using the orthogonal approximation are derived in Appendix A. They 
are extensively used in Sec. VI in the calculation of the critical indices using the orthogonal 
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approximation. 

Section VII describes the calculation of the critical exponents for the isotropic cases 
utilizing the orthogonal approximation. The corresponding loop integrals are computed in 
Appendix B. 

Section VIII presents the exact computation of the Feynman integrals in the isotropic 
behaviors. The reason for this explicit computations is that the four-point graphs for arbi- 
trary n is quite difficult to get in a closed form explicitly. By fixing n we can compute their 
contribution and find a recursion formula for arbitrary n. In addition, the two-point vertex 
part graphs are also computed up to three-loop order. They are shown to be simpler that 
the four-point contributions. 

Section IX is an exposition of the exact results obtained for generic isotropic critical 
exponents for arbitrary n. They are determined diagrammatically using the results of Section 
VIII. They are shown to be identical to the exponents previously evaluated using the massless 
formalism. 

We present the discussion of our results along with the conclusions and further possible 
applications of the present method in Sec.X. 



II. NORMALIZATION CONDITIONS FOR THE MASSIVE THEORIES 



The functional integral representation of the CECI model was first introduced in Ref. 35|. 
It corresponds to a theory containing higher derivatives. The larger the number of 
neighbors coupled via competing interactions, the higher is the power of the derivative 
terms in the bare Lagrangian. In the anisotropic behaviors for generic competing systems, 
there are many simultaneous types of higher derivative terms. The original bare Lagrangian 
density reads 

L = ^IV(.-sL.™.)<^ol' + ETlvL0or (1) 

n=2 



L ^ L—1 n— 1 ^ 

n=2 n=3 n'=2 

The first and third summations in the above Lagrangian correspond to the effect of the 
competition occurring in the system. The parameters cr„ and r„„/ are introduced to guaran- 
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tee that all terms have the same canonical dimension. At the Lifshitz critical region, the fixed 
values of the exchange interactions are transliterated into the conditions 6on = Tmi' = in 
the above Lagrangian. This simplification permits the decoupling of each subspace in Feyn- 
man loop integrals, but produces the apparent complication that all the higher derivative 
terms become relevant in the free massive propagator. Notice that the Lifshitz point is 
characterized by //q = when the temperature is exactly at the Lifshitz value T = Tl, with 

Son = Tnn' = 0. 

It is worthy to separate the momentum subspaces in the form = = 
(p, k2, Icl) where = p is a vector along mi spatial directions connected to the first 
(non) competing subspace, p(2) = k2 where k2 is a vector along m2 space directions asso- 
ciated to the second neighbor competing subspace, etc., = where kL is a vector 
along rriL directions associated to the L-th neighbor competing subspace. The variation of 
Kn in the renormalized theories whose starting point is the same bare theory are induced 
by the existence of independent correlation lengths These L independent flows in the 
momenta can be implemented through L independent renormalization group equations for 
each m„- dimensional subspace. Nevertheless, this construction is consistent since the appar- 
ent overcounting producing L independent coupling constants defined in each spatial subset 



can be overcome, for all of them f 
As discussed in detail in Ref. 



ow to the same infrared nontrivial fixed point. 



361] in the massless theory, the conditions which define 



the critical region can be used to perform a dimensional redefinition of momentum scales 
along each type of competing subspace as follows. Let [p] = M be the mass dimension 
of the quadratic momenta corresponding to the noncompeting subspace mi. We can get 
rid of the parameter (T„ appearing in front of each higher derivative term characterizing 
each competing subspace by redefining the associated momenta scales through the relations 
[k(n)] = Mn. This effectively disentangles each momenta subspace such that they can be 
treated independently. In the anisotropic criticalities there are L subspaces, each of them 
are m^^-dimensional. The masslessness of the theory at the Lifshitz temperature requires 
that the renormalized IPI vertex parts (the basic objects in this framework related to the 
thermodynamical potentials of the critical system) must be computed at nonzero external 
momenta in order to avoid infrared divergences. For instance, if we consider the vertex 
functions along the L-th subspace we use normalization conditions by choosing the symmetry 
point at nonzero external momenta kl along this subspace, whereas all the other momenta 
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scales perpendicular to the mi directions are set to zero. Therefore, the nonzero momenta 
used to renormalize the theory in each individual subspace can be viewed as a label in 
parameter space defining L independent sets of vertex functions. On the other hand, there 
is no need to consider more than one subspace for isotropic behaviors. 

Now, let us take a look at the simplest CECl model, namely, we take m2 = = 1 with 
rUn = for n > 3. The competition are located at the y and z axis. The phase diagram 



(see Fig.l of Ref. 



36|]) can be described by T,pz = J2z/ Jiz,Piy = J2y/ Jiy and p2y = Jsy/ 



ly 



We can represent the phase diagram in a very simple form through several two-dimensional 
projections by fixing some parameters and varying only two of them. For instance, the 
diagram (T, p^) is associated to the second character behavior provided piy and p2y are kept 
fixed. On the other hand, the purely third character behavior can be expressed in terms 
of a three-dimensional phase diagram with axes {T,piy,p2y) so long as pz is fixed. But 
this can be further simplified by considering only the two-dimensional projection of this 
phase diagram as suggested above if we take instead the variation of {piy,p2y) with (T,pz) 
at constant values. Now the superposition of these two two-dimensional phase diagrams 
yields a point of intersection for particular values of the temperature which is identified with 
the generic third character Lifshitz point. Notice that the ferromagnetic phase and the two 
modulated phases named Helical2 and Helical^ (see Fig. 3 of Ref. |3Q]) encounter themselves 
at the uniaxial generic third character Lifshitz point. There are two first order lines: one of 
them separates the ordered- if eZica/2 regions whereas the other splits the Helical2 — Helical^ 
phases. 

This situation can be generalized for generic higher character Lifshitz points if we split the 
corresponding multiparameter(multidimensional) phase diagram in two-dimensional slices 
and by superposing them together in a single two-dimensional diagram. There are now L 
modulated phases meeting at the L-th generic higher character Lifshitz point. We emphasize 
that each competing subspace is defined by its own independent correlation length, i.e., ^1 
for the subspace with only ferromagnetic exchange forces coupling first neighbors, C,l for 
the subspace defined by alternate signs in the exchange forces up to L-th neighbors. 

Roughly speaking, we can associate the (inverse of the) mass to the correlation length. 
The anisotropic generic higher character universality classes require L independent correla- 
tion lengths. In close analogy to what has been done for the m-axial Lifshitz critical behavior 
(actually a second character Lifshitz behavior) we can reexpress the bare Lagrangian density 
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(1) in terms of L independent bare masses as follows. 

It is obvious that we can attain the L-th generic higher character Lifshitz point in the 
phase diagram outlined above by varying the "mass" coming from the noncompeting sub- 
space. This means that we approach the L-th generic higher character Lifshitz point from 
the ferromagnetic phase. Simple inspection of the phase diagram indicates that we can reach 
this multicritical point coming from any of the several modulated phases. Then, it is possi- 
ble to introduce L independent bare masses such that they generate the L renormalization 
group flows in parameter space which are compatible with the L correlation lengths present 
in these criticalities. Let us describe the introduction of the many independent mass scales 
beyond these simple phenomenological considerations. 

Siegel's method of dimensional reduction suggests in a simple way how to introduce mass 
in ordinary quantum field theories with quadratic derivatives in the Lagrangian density 



The basic steps for scalar fields are: extend the range of the momenta indices and call the 
extra index choose the momentum component associated to this direction equal to the 
mass p_i = /io and introduce factors of i to re-establish reality di = ip-i = ifiQ. The operator 

in the higher dimensional space (including the extra index results in a massive 

operator + /Xq in spacetime dimensions. In the Lagrangian (1) the metric is Euclidean 
but the conditions (5on = Tnn' = can be used to perform the dimensional redefinition of the 
momenta characterizing each type of competing axes. The dimensional redefinitions turn out 
to implement independent dilatation invariance along the mi = d — Yln=2 ™n noncompeting 
directions, m2 subspace, and so on, up to rriL space directions. The typical momentum 
combinations which appear in the inverse free propagator is of the form pf + Yln=2i^n)"'- 

The extension of the massive method to include arbitrary types of competing axes can be 
understood from the analysis of m-axial critical behavior employing massive fields. The in- 
troduction of distinct masses using Siegel's recipe translates itself in the following conditions: 
extend the range of the vector indices in the n subspace to the "extra direction" " — 1"(^), 
identify the conjugate momentum component with the mass in that subspace P-i(n) = /^on 
and use factors of i to restore reality (9_i(„) = ip-i(n) = ifJ^on- If we apply this reduction 
to the n = 1 subspace, the operator pf + Ylt=2(^n)"' the higher dimensional space turns 
out to become pi + J2n=2i^n)"' + /^oi ^ space dimensions. When we utilize the same 
procedure to the subspace n = 2 with /ci(2) = /io2, the simplest situation occurs for uniaxial 
case 1712 = 1- The combination in the inverse free propagator becomes Pi + + Ylt=3(^n)^- 
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If we define the internal product in the higher dimensional space with index " — 1(2)" by 
^2 + ^-1(2)' dimensional reduction yields + + J2n=3i^lT ~^ Pi + ^2 + J2n=3i^n)"' + /^02- 
Sticking to simplicity, we can extend this procedure to the n' subspace in the uniaxial case 
rrin' = 1- Choosing the internal product in the higher dimensional space including the 
index —l{n') as k'^?' + the combination pj + k'^f + J2n{^n')=2(^n)"' is reduced to 

Pi + ^n"' + X]n(^n')=2(^n)" + lAn' ■ Thcsc arguments suggest that it is possible to choose the 
mass defining each competing subspace with the same canonical dimension of the momenta 
along those directions. 

The resulting strategy leads us to define the bare masses in the bare Lagrangian density 
with different powers, depending on the chosen subspace we work with. This makes explicit 
reference to the fact that those distinct subspaces are inequivalent. Our experience handling 
labels in the massless theory suggests that when dealing with the renormalized theory, the 
bare mass in each subspace naturally defines a renormalized mass and coupling constant 
characterizing that subspace. We can then implement L independent bare masses along 
with L independent bare coupling constants in the bare Lagrangian density. Bearing in 
mind these considerations we write the original bare Lagrangian in the form 

^ = ^IV(.-j:^.,™„)0oP + X:YlvL</'or (2) 

n=2 

L ^ L—l n—l ^ 

n=2 n=3 n'=2 

We shall focus our attention hereafter in the Lifshitz critical region where 5on = Tnn' = 
in Eq.(2). The independent bare mass in each competing subspace naturally prevents the 
appearance of infrared divergences in calculating Feynman integrals for the associated IP/ 
vertex functions. Thus, the label n can be used to define independent renormalized IP/ 
vertex functions which have the following property: those with arbitrary external momenta 
along the n-th space directions have a nonvanishing bare mass /^on and coupling constant 
Aora. However, they have vanishing external momenta {pn' = 0), bare mass (/ion' = 0) and 
coupling constant (Aon' = 0) along all directions elsewhere (if n' 7^ n). Moreover, the 
zero mass limits of the vertex parts are well-defined and reduce to the cases previously 
investigated in Refs.|35l. |36[|. 
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One important ingredient to complete the description is the existence of L independent 
cutoffs, which are responsible for the independent variations in the mass parameters of the 
anisotropic cases. Each cutoff has the same canonical dimension as mass and momenta 
characterizing the competition subspace under consideration. There is a similarity and a 
difference when we compare the massless and the massive approaches. Unlike the massless 
case, the massive theories do not flow to fixed points where they become scale invariant: 
scale invariance is only achieved if the coupling constants in each subspace are set exactly 
at the eigenvalue conditions Unoo, i-e., at the non-attractive fixed points. On the other 
hand, the ultraviolet fixed point value of the coupling constants m„oo are independent of n, 
a feature already encountered in our analysis of the massless theory in the infrared regime. 
The isotropic behaviors can be described by the Lagrangian (2) with some modifications: 
the first term does not appear in the isotropic situation d = rrin, n = 1, 2, .., L and there is 
solely one kind of bare (renormalized) mass. 

We begin by describing the anisotropic renormalization conditions. Except for minor 



modifications, we shall follow the conventions adopted in Refs. |31l . I36l ] . Let the n-th subspace 
(n = 1, L) be defined by a nonvanishing value of bare mass and coupling constant /ion 7^ 
0, Aon 7^ with ^Qn' = 0, Aon' = for n' 7^ ra. The nonvanishing bare parameters induce 



renormalized parameters /i^ and The slight change of notation with respect to Ref. [31| 
in defining the renormalized mass is performed in order to prevent the confusion with the 
number of space directions m„ of the n-th competing subspace. The renormalized IP/ 
vertex parts living in the m^i-dimensional subspace are defined by the following normalization 
conditions: 

r|^„)(0,/i„,^/n)=/i^", (3a) 
d^n k)=o = l' (3b) 

rK(„)(0,/in,^n) = ^n, (3c) 

rg;^5(0,0;0,/i„,(7O = 1. (3d) 

In analogy to the massless case, we can fix the renormalized mass scale in each competing 
susbspace by choosing /i^" = 1. 

The isotropic situations can be described directly by using the above normalization con- 
ditions for the renormalized vertex parts where only one type of competing axes are required 
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in that case. Although the scahng analysis is a bit different in the isotropic and anisotropic 
cases, they share the same normalization conditions Eq.(3). 



III. ONE-LOOP RENORMALIZABILITY AT THE CRITICAL DIMENSION 

In this section we shall investigate the divergence structure of graphs corresponding to 
primitively divergent bare IPI vertex parts. In general, the divergences can be expressed 
in terms of regularized expressions which retain their infinite values provided the regulators 
used take appropriate hmits. We shall show that these divergences can be expressed rather 
simply in terms of independent cutoffs in the anisotropic cases, when there are many in- 
dependent competing subspaces appearing simultaneously in the problem. Otherwise, the 
isotropic cases only require one type of cutoff. We discuss the renormalization of all ver- 
tex parts at one-loop level. Many results of this section are going to be useful in proving 
renormalizability at arbitrary loop order. We shall restrict our attention throughout only 
to vertex parts which can be renormalized multiplicatively. 



A. Anisotropic Sector 

Consider the noncompeting subspace corresponding to the label n = 1. The vertex parts 
associated to it possess arbitrary external momenta Pi = p along the mi = d — X]n=2 
space directions, with nonvanishing bare mass ^oi and coupling constant Aqi. The bare 
primitive divergent vertex parts are rj^(i) , ^bIi) ^^(i) ' '^^^^^ perturbative expansions 
up to one-loop order can be written as: 



{l)\Pi) - ^01 ^( 



[EL2ikUr + i<l + Pi + + (Et2(^?n))" + ?'+/^0l) 



+ {Pi Ps) + (P2 Pa)), (5) 

^ ' 2 ./ (12 ^r^ I ('^ I ^ ^2 I ,,2 1 / V^i (1.2 ^r^ i „2 i ,,2 \ 



[T.n=2K)Y + +P3)^ + l^ll] {T.n=2{kl)Y + + l4.) 
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It is easy to see that for m2 = ... = mi = the above integrals reduce to the ordinary 
noncompeting ones and there is no subintegral involving competing momenta. However, 
we shall keep the values of m„ unspecified and compute the integrals at the upper critical 



dimension 4 = 4 + Y.n=2[^^^ 

We begin with the computation of the integral with a single propagator appearing 
in Eq.(4). We may employ the Schwinger's trick to write the propagator in terms of a 
parametric integral, namely 

-j^ poo L 

^- ^= / da exp[-aiY,ikfn)r + q"" + fill)]. (7) 



Making use of the identity [27 1 

j dxi...dxm„exp{-a{xl + ... + x^^)") = ^r(^)a^^„„, (8) 
we can solve the integrals over the momenta. The resulting expression for Ji(i) is given by 

= 2VEL.^)r(2-El7)(nL2 "'\ ''' ) / exp[-af^l,]a-\ (9) 

n=2 

We introduce the cutoff Ai in order to express mathematically the ultraviolet divergence 
implicit in the parametric integral in terms of it. The divergence comes from the region for 
small values of a. We then regularize the integral by suppressing a domain (0,A^^) in the 
small a region of integration. It is very simple to introduce in the integrand a function of 
a and Ai /a^ (a) which vanishes for a < A^'^ (whose derivatives vanish in the limit a — 0) 
and is identically one for a > Aj"^. The Heaviside function is efficient to produce this effect 
and we choose /ai(«) = ^(a — A|f^). Integrating by parts twice and letting the cutoff go to 
infinity, the divergence can be expressed in terms of the cutoff as 

The overall angular factor is different from the pure 0^ field theory and appears in the same 
way in both cases: it shows up whenever a loop integral is performed. Proceeding in the 
standard way, it can be absorbed in a redefinition of the coupling constant. Looking at the 
singularity structure of we find that its dependence on Ai is exactly the same as its 
counterpart describing ordinary critical behavior has in terms of the cutoff, say A, at the 



critical dimension d = 4 
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The integral contributing to both T^^^^ and T^^^^j denoted by /2(i), namely 

can be performed in a similar fashion. Introduce two Schwinger parameters and integrate 
over the momenta. We are left with two parametric integrals over ai and a2. Integrate first 
over «i by defining a new variable a' = ai + a2. The integration limits of the integral over 
a' turn out to be {a2, oo). This integration produces no divergence, but the integral to be 
done over a2 results in the ultraviolet divergence for small values of a2- We regularize this 
integral exactly as before by introducing the cutoff function /ai(q;2) = 0{a2 — A]^^) inside 
the integrand. Expanding the integrand in powers of the external momenta, we find out 
that the divergence is present only in the momentum independent term and taking the limit 
Ai — > cxD, we obtain the following divergent result 

L 



Owing to the similarity with the m-axial Lifshitz situation, when there is only one type 
of competing axes, we are going to treat all competing subspaces at once. We select only 
the competing axes in the m„-dimensional subspace where competing interactions couple 
up to n neighbors. In that case, we start with /io„' = /Uon^nn'j '^on' = ^on^nn' and external 
momenta Pi{n') = ^i(n)<^nn'- The vertex functions which have primitive divergences can be 
expanded up to one-loop order as 

rS(^H) = (Cr + .^^+^ / f f-^fe-^-^H , (13) 

2 7 (j:L2ikur + + i^in) 



^7 } ~Y + (K(n) ~^ K(n)) + (K(n) ~^ K(n)))^ (14) 

{zL2ikur + + ^'c;^) 

(2,1), w w Ao„ f d'^-E^=.m„gn^^,d"^-fc(.) 



p(2,l)/,/ ,/ ,/ N _ -, Aon f 
^ (n) I'^lW' '^2(n); %(n)J — 2 / 



2 J [EfnYn)n'=2(^?n'))"' + ((^W+^3(n))')" + ?' + /^in] 

X7 ^- (15) 

(Et2(^?n))" + ?^+/^gn) 
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In order to perform these integrals, we recall that the bare mass parameter has ^ of 
the canonical dimension of the momenta along the noncompeting (quadratic) subspace. We 
then choose the cutoff A„ associated to this subspace with the same canonical dimension 
as /ion- When using the Schwingcr parameters to perform the integrals over the momenta 
variables, we encounter the divergences implicit over the parametric integrals in the region 
of small values of a. We cutoff the parametric integrals introducing in the integrand the 
regularization function /a„(q;) = 6{a — A~^"), which restricts the integration domain to the 
A~^" < q; < oo. Denote the integrals contributing to the two-point function at one- loop by 
Ii{n) whereas J2(„) represent the contributions to v'^^^-^ (and rl^'^^"*), respectively. Taking the 
limit A„ — >■ CX3 we get to the following expressions for these objects {n = 2, L) 

1 m ^ r(m:iL\ A2n A 2n 

1„ Am„,._. SmSi^)., .A?" 



Comparing the above equations with their counterparts in the noncompeting (quadratic) 
subspace, it is obvious that we can now unify all subspaces by taking n = 1, ...,L in the 
above formulae. This unification will be useful in our description of the renormalization of 
masses and coupling constants at one-loop order. The fact of the matter is that when we 
use the normalization conditions at zero external momenta in each subspace with a set of 
finite renormalized parameters (//„, Qn) starting from infinite bare quantities {/lon, Aon), we 
can express the former in terms of the latter via the following equations: 

/^r = /^On + ^/l(n)(0), (18) 
yn = Ao„-^/2(n)(0). (19) 

If we go the other way around by writing the bare parameters in terms of the renormalized 

ones and getting rid of higher order corrections in the renormalized coupling constants, these 
simultaneous operations render the original bare vertex parts T^^^ and T^^^ finite. In fact, 
the corresponding finite vertex parts are given by {n = 1, L) 

rS(P(n)) = (P(n))^'^ + /^r, (20) 
^'S)^Pi{n)) ^9n- y(-^2(n)(Pl(n) + P2(n)) + -^2(n) (Pl(n) +P3(n)) + -^2(n)(P2(n) +P3(n)) " 3/2(n)(0)). 

(21) 
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The bare vertex parts T^^^j with N > A have skeleton expansions. There are L independent 
sets of skeleton expansions, one set for each competing subspace. Thus, their diagrammatic 
expansions result in finite expressions at two-loop level when they are written in terms of 
the renormalizcd mass(cs) and coupling constant (s) at one-loop order. 

In order to study insertion of composite operators and their renormalization, we analyze 

(2 1) 

the bare vertex F^^'^ '. They are not automatically finite by the reparametrizations turning 
the bare mass and coupling constant into those renormalized amounts. We can define the 
renormalized (finite) vertex part T^^^^^ as 

An). (22) 

The normalization conditions (3d) require that, at this order in the loop expansion, the 
normalization functions are given by 

V(n) = l + f^2(n)(0). (23) 

Replacing this expression into the definition of T]^^^^ turns out to make this vertex function 
finite, which may be written as 

r§n) ' P2(n) ; P3(n) ; Qn, A^n) = 1 - y {hin) {P3{n)) " h{n) (0)) . (24) 

In most applications we shall describe insertion at zero external momenta, i.e., we refer to the 
vertex function r^(„)(P(n), —P(n)', 0; Qn, l^n)- In this manner, we can address the multiplicative 
renormalizability and obtain independent fiows in the parameter spaces of vertex parts 
including arbitrary types of composite operators, irrespective of the competing subspace 
under scrutiny. 

The remaining multiplicatively renormalized vertex parts including composite operators 
'^^R(n) with (A^, L) > (2, 1) have skeleton expansions. Consequently, they are finite at two- 
loop level whenever we use gn, f^n and T'^^l^^{p^n),P2(ny,P3(ny, 9n, l^n) (or Z02(„)) inside their 
one-loop subgraphs, due to the following high momentum pattern of the primitively divergent 
vertex functions 

\^Rln)(PnP{n))\ < Pn" X POWer of lupn, 
\^Rln)iPr^Pi{n))\ ^ POWCT of lupn, 
\^Rln){PnPl{n),PnPl{nyPnq(n))\< POWCT of lupn, 
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at every finite order in the limit p„ — > oo. This large momentum behavior is not going 
to be derived here, but we shall assume that it is valid henceforth. In other words, the 
Born values of the various renormalized vertex functions are modified only by powers of 
logarithms. These modifications are caused by the interactions in every perturbative order at 
the loop expansion, and we can use a power counting reasoning to approach the multiplicative 
renormalizability at the critical dimension to all loop orders. This (nonperturbative) proof 
of multiplicative renormalizability corresponding to the Callan-Symanzik-Lifshitz massive 
method shall be analyzed after our construction of similar renormalization background for 
isotropic behaviors at one-loop order, to which we turn our attention next. 



B. Isotropic Sector 

Recall that for arbitrary higher character isotropic Lifshitz critical behaviors there is only 
one type of competing subspace. The critical dimension of the isotropic critical behaviors 
when there are n neighbors along each space dimension interacting via alternate exchange 
forces is d — m — An, with n — 1,...,L. The one- loop required vertex parts which are 
primitively divergent are given by 

\ r Jin h.1 

rSw = (w-r..g;.^ . (25) 



Let and /2(n) be the one-loop integrals involved in the calculation of F^J^ and r|J^, 
respectively. The computation of the integrals are entirely analogous to the nth competing 
subspace of the anisotropic sector: the cutoff A„ has the same canonical dimension of //on 
and the removal of small values of the parametric integrals in the Schwinger parameters are 
implemented via Heaviside's function. The only difference with respect to the anisotropic 
case study in the nth competing subspace is that the angular factor is different. We then 
obtain 

-I \2n A 

^Hn){k) = - ^4n f^l:i-p,-H-^)), (28) 
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hin){k) = ^ 54n ln&. (29) 

Renormalization of the mass and coupling constant can be accomplished through the uti- 
lization of the normalization conditions (3a) and (3c). When expressed as functions of the 
bare amounts (which have infinite values) they take the form 

/^f = /^On + ^/l(n)(0), (30) 
^n = Ao„-^/2(„)(0). (31) 

We could have written the above expressions in the opposite direction. In that case, either 
the bare functions F^f^^ (ki){iJ,n, Qn) can be rendered convergent at one-loop order provided 
N < A, or they are skeleton expansions. To conclude, whenever we use the normalization 
condition (3d) along with the appropriate version of Eq. (22), which amounts to define the 
normalization function Z02(„), they produce a convergent expression for the renormalized 
vertex r^^^^^(A;i, k^). Consequently, this leads to Z^2(^n) = 1 + ^-^2(n)(0). 

Next, let us tackle the issue of multiplicative renormalization in order to extend the one- 
loop discussion just explained in the present section to arbitrary loop order employing the 
Callan-Symanzik-Lifshitz equations. 



IV. THE CALLAN-SYMANZIK-LIFSHITZ EQUATIONS AT THE CRITICAL DI- 
MENSIONS 

Let us establish the inductive proof of multiplicative renormalizability and prove the 
existence of the Callan-Symanzik-Lifshitz equation in the presence of various simultaneous 
types of competing axes as well as when solely one type of competing axes take place in the 
generalized Lifshitz critical behaviors. 

Although we shall not make explicit reference, we shall be concerned with vertex parts 
which arc regularized using cutoffs as described in the previous section. It is implicitly 
assumed that this is the case in the remainder of the present section. Our discussion here is 
restricted to the critical dimension of the theory. 
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A. Anisotropic behaviors 



We are going to study the physical system at its upper critical dimension dc — 4: + 
X]n=2[ ^^n^^ ]^"- multiplicatively renormalized vertex parts including composite opera- 
tors are defined by 

^Bln)\Pi{ri),Qi{n),gn, l^n) = Z^^^^Z^2(^n-)T^'^'^\pi(^n), Qi{n), ^On, I^On, ^n) ■ (32) 

We emphasize that vertex function with (N, L) = (0, 2) which arc additively renormalized 
are precluded from our analysis. In the above expression, (i = 1, ...,N) are the external 
momenta associated to the N external legs of (p operators, (5i(n) — 1? the external 

momenta corresponding to the L insertions of 0^ composite operators and the independent 
cutoffs A„ are required to depict the regularization process in each inequivalent subspace. 
In the present section, we shall suppose that every bare diagram at any loop order can be 
implicitly regularized with the appropriate cutoff. 

When we apply the operation to the bare vertex part T''^^f\pi{n), Qi{n), ^On, l^-On, An) 
with fixed Aon and A„, we obtain a zero momentum insertion of the operator in the vertex 

P^^* ^[n)^\Pi{n),Qi{n),>^Qn,I^On,K), namely r\^'^j^'^^\pi{n),Qi{n),0,XQn,I^Qn,K)- Then, WC 

have 



Q (ATX,) (7VL+1) 
rj 2n^{n) Qi(n); -^On, /^On; A.„) = r. ^' (Pi(n) ) Qj(n) ; 0, Aon, A'-On; A^) . (33) 

C/ig,„ 



Eq.(32) can be employed to transform the bare quantities into the renormalized amounts. 
Writing the bare mass //on representing the m„-dimensional subspace in terms of the renor- 
malized mass and coupling constant (//on — l^'Onil^'n, 9n)) a-nd applying the chain rule, the 
renormalized vertex functions obey the equation 

+ - L^)^'n(S^P^in), 9n, ^^n) = (34) 

p(Ar,L+l)/ ^ p. X 

R{n) \Pi{n) ) Vi(n) )'^) 9n) IJ^n) i 

Where 2np„ = ||V(n), ^ = V(n)|%, ^ = Vh^, = Vh^^- K 
we define the functions /3„(= ^) = 7^(n)(= ^) = /^n^^^ and 7</.^(n)(= -f^) = 

~l^n—d^^, together with the multiplication of the last equation by ^ yield the result 

S/T y ^'/'C") + -^7</.2(n) ) r^;^^^ (Pi(n) , , gn, /^n) = (35) 



0^..2n ^Mo" /7-I p(iV,L+l)/ ^ -n ^ 



21 



In order to express last equation purely in terms of renormalized objects we make use of 
the normalization conditions r^(„)(0) = /i^", ri^(n)(0) = 1 the particular case {N,L) = 
(2, 0). This results in the Callan-Symanzik-Lifshitz equations (CSLE) for anisotropic Lifshitz 
points of generic competing systems: 



^^"''dfT ^"9^ y^'/'W + ^'y4>^{n))'^R(n) {Pi{n), Qi{n), 9n, ^^n) = (36) 



Some comments are in order [2l|. Initially, apply a derivative of the CSLE at zero 
external momenta for (A^, L) = (2, 0) and use the normalization conditions of section II. 
This means 



d 

7<^(n) = /xf (2n - 7<^(n))70^rg^5(p(„), 0, gn, /in)|p2« =o- (37) 



(2 1) 

At one-loop order we already encountered from Eq.(24) that r]j^^j(p(„), — p(„); 0, /i„) = 
1 + 0{g^). It is easy to see that 'y^{n) begins at 0{g'^). Next, substitute {N,L) = (4,0) in 
the CSLE. Taking advantage of the normalization conditions leads to 

f3n - 270(n)^n = /in"(2n - 70(n))r2(i) (0, 0, 0, 0; 0, fin, ^n) • (38) 

The first contribution to r^^^|(0, 0, 0, 0; 0, (?„) is 0{g^). From last equation we discover 
that (3n is 0{g'^). Consequently, the operation /3n-^ is 0{gn), differently from the operation 
which is 0{g^). Now, taking N = 2 and L = 1 into the CSLE along with the 
normalization condition r]:j^Jj(0, 0; 0, /i„) = 1, we find 

- 7</.{n) + 702{n) = /^n"(2?^ - 7</.(n))r2(^5 0; 0, 0, fin, gn) ■ (39) 

(2 2) 

As 

^R{n) starts at 0{gn), we find out directly that 7(/,2(„) = 0{gn)- We have at hand the 
basic requirements to pursuing the inductive proof of multiplicative renormalizability at all 
orders of perturbation theory. We turn now our attention to this subject. 

We begin the inductive proof of multiplicative renormalizability with the hypothesis that 
the renormalized vertices defined by Eq.(32) have been transformed into convergent expres- 
sions up to the Lth-\oop order at fixed /i„ and gn in the limit A„ — > oo. Thus, we commence 
with the claim that in the infinite cutoff limit fI?) ^ Fp? ^ and T%^\ are finite at order gff', 

H(n)^ K(n) H(n) ' 

5'^"''^ and g^, respectively. We suppress the arguments of the renormalized vertex functions 
in order to simplify the succeeding discussion. 
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The preliminary observations considered above at one-loop can be generalized to L-loop 
level. Using Eq.(37) at this order in the coupling constant gi„, we find out that 70(n) is finite 
at O(g^). Prom Eq.(38) at (L + l)-loop order, T^^^^s^ in the right-hand side (rhs) is already 
finite at (L + l)-loop order for it has a skeleton expansion. On the other hand, at (L+l)-loop 
order this vertex part is of order gn~^^- In other words, the combination — 2j(f,(^n)gn) is 
convergent at 0{gl^~^^). As 7</,(n) is finite at 0{g^), this implies that /?„ is finite at 0{g^~^^). 

We can write the CSLE alternatively as 

(/^n^^) = + f 7.(n) - M7,.(„))rg[;^^ + (2n - 7.H)rSJ[^^"^^ ■ (40) 

Our aim is to prove that we can find a finite result for the rhs of this equation at (L-|- l)-loop 
order. Note that T^^^-* in the rhs is only needed at Lth loop order, since its coefficients 
Pn-^ and 702(„) are at least of order gn- In the last piece of the rhs, either T^^^^^^ has a 
skeleton expansion (convergent at (L + l)-loop order) or the CSLE must be iterated. 

Let us fix our attention in the case = 4, M = in Eq.(40). Both terms in the rhs 
are finite at (L -|- l)-loop order. We then conclude that {t^-n—gj^) is also convergent at this 
order. 

The proof that g^^"' ) is finite at (L-|- l)-loop order as well can be achieved by consid- 
ering the perturbative integration of L^^^^^. At the critical dimension L^^^^^ is dimensionless; 
it gets unchanged under independent dilatations in their dimensionful parameters: 

(Pn, /^n, A„, gn) ^ (PnPn, Pn"^n, Pn^n, 9n) ■ (41) 

The choice Pn — together with the dilatation invariance stated above can be written as 



r2(n)(Pn, /^n, K: Qn) = ^f(n)(^^ 7^' 9n)- (42) 

P^n Pn 

Taking = 4, M = in Eq.(40) for arbitrary external momenta at order g^'^'^i it can be 
reexpressed in the following manner in terms of a running variable //^ 



/ ,^^R{n),Pn A„ M)/Pn A„ 

^^'r^^r;^^^^^^9n))\L+2^fl^~^{ — ,—,gn))\L+2. (43) 



First, write the running mass as a function of a dimensionless variable /i'^ — The 
integration of /i'^ in the interval [oo, amounts to integrate over a in the region (0, 1). We 
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then employ the normahzation conditions (3c) at the hmit //[j = oo as boundary conditions. 
We encounter the following result 

rT-(4) /Pn xn, da . (A) Pn A„ , . 

[ry„)( — ,— ,^n)]|L+2 = ^n- / —[/(^ (« — ,« — , ^„)] |l+2. (44) 

From what we have been discussing, [/^''^j]|L+2 can be written entirely in terms of lower 
order amounts which, by assumption, all possess a finite value at A„ — > oo. Since [/(^j]|L+2 is 
analytic for small momenta, the limit of zero external momenta can be taken safely. Indeed, 
this does not introduce any difficulty into the integral over a. Hence, we are led to 



[rR(„)( — ,00,5n)]|L+2 = ^n- / — /(;)'(« — , OO, ^„) |l+2, 
t^n Jo " H'n 



(45) 



which demonstrates that the renormalized vertex F^^^^^ can be successfully related to a finite 
integral over lower order renormalized vertex functions. 

Considering Eq.(40) for the values A?" = 2, M = 1 at order g^'^^ yields 



(46) 

The first term in the rhs is a combination of /3n at 0{g^^^) and F^^^^^^^ at 0{g^), for Pn-^ 
is 0{gn)- The two contributions coming from this term are finite at this loop order in the 
above equation. The last term has to do with F^^^^^ which has a skeleton expansion. Since 
the couphng constants and masses are actually finite at 0{g^), this implies that the last 
term is automatically finite at 0{g^'^^). 

The second term in Eq.(46) has two contributions: T^^ll^^ at 0{g^) which is hypothetically 
finite (7^2 („) starts at 0{gn)) as well as (j^n) - 7<^2(n)) at 0{g^+^). The proof that (7^(„) - 
702(-„)) is finite at 0{g^^^^) still has to be demonstrated. The calculation of this expression at 
zero momentum is identical to Eq.(39) computed at (L + l)-loop order. The rhs of Eq.(39) 
involves the skeleton F^^^^^ at this loop order, which is finite. This shows that (7</,(n) — 702(n)) 
is finite at 0(51^+^). Consequently, this concludes the proof that the rhs of Eq.(46) is finite. 

Integration of the left hand side of Eq. (46) for arbitrary external momenta gives the result 

^r(2,i) 

(/..^)U.. = 4r ■ (47) 
(2 I) 

Owing to the dimensionlessness of F]^^^^^ at the critical dimension, we can perform a scaling 
transformation following similar steps to what was done for the vertex F^^^^^. By requiring 
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analyticity of f^^-^ ' for small p„ and using the normalization conditions at zero external 
momenta, the proof that T^^^^^ is manifestly finite at 0{g^~^^) follows at once. Using Eq.(37), 
we learn that 70(„) is finite at 0{g^~^^) as well. 

The conclusion of the inductive proof of multiplicative renormalizability can be achieved 
by turning our attention to the case = 2, M = at order gf^+^ in (40), which reads 
5r(2) ^ 

^^"^^'"^^ = + [7^(n)rg(:5]|.+x + [/.r(2n - 7^H)rg:5]U+i . (48) 

Now, T^^jn) required solely at 0{g^J (finite by hypothesis) since /3„ starts at 0{gn) and we 
have already proven that and 70(n) arc finite at 0{g^^^). This is sufficient to prove that 
the first two terms of the rhs in the above equation are manifestly finite. For the last term, 
we have just proved that both contributions inside it are finite at 0{gl^~^^). Since the rhs 
has a finite limit, we proceed to integrating the CSLE for this vertex function with arbitrary 
momenta. Perform the redefinition 

rg(l)(Pn) = rg„)(p„)-pf-A.r. (49) 

The redefined vertex has mass dimension and only deviates from r^^^^-|(p„) by higher 
powers of pf^ and fif^, which turn out to be cancelled by negative powers of the cutoffs. 
Normalization conditions require that it vanishes with {p^Y ^^^^ small \pn\- Thus, collecting 
all the information contained in our previous discussions, T]^^^^^ (p„) satisfies an equation of 
the type 

~ (2) 

with = 0((p^")) for small Introduce the change of variables ji'^ — Integrate over 
the variable ^'^ in the interval (oo,//„). Using the normalization condition at zero external 
momenta as a boundary condition to the solution, i.e., r^^7j^(0, 0, (7^) = 0, and taking the 
limit of infinite cutoffs, we get to 

We are left with the task of proving that the integral is finite. But this is straightforward 
from the behavior of the integrand for small it goes like /(f)(ct^) CC) 9n} — 0((ap„)^")^), 
which is free of singularities in the lower integration limit a — > 0. 

We have thus succeeded in demonstrating the multiplicative renormalizability by the 
inductive method as well as the existence of the Callan-Symanzik-Lifshitz Equations (36). 
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B. Isotropic behaviors 



The isotropic situation has a simple parallel with the nth competition subspace of the 
anisotropic case as explicated before. We just have to keep in mind that the critical dimen- 
sion is different dc = and there is only one type of competition along all space directions. 
Then it follows in a straightforward manner that the isotropic vertex parts satisfy exactly 
the Callan-Symanzik-Lifshitz equation associated to the nth competing subspace of the 
anisotropic behaviors discussed above, namely 

^ ^"^dg y^'^^") ^'y<t>^{n))'^^R(!^\pi{n), Qi{n),9n, /^n) = (52) 

In addition, the arguments presented above in the inductive proof of renormalizability for 
the nth subspace in the anisotropic cases can also be used to prove inductively the finiteness 
of all vertex parts of the isotropic cases. The situation is identical to the proof furnished 
in the particular isotropic m-axial case (n = 2) Jsil. By the same token, the existence of 
the Callan-Symanzik-Lifshitz equations (52) for arbitrary isotropic critical behaviors follows 
directly. This concludes the formal proof of multiplicative renormalizability of arbitrary 
competing systems of the Lifshitz type to all orders in perturbation theory. 



V. THE CALLAN-SYMANZIK-LIFSHITZ EQUATIONS AT d = 4 - 

We will restrict our considerations only to vertex parts including composite operators 
which can be renormalized multiplicatively. The renormalized IPI vertex parts are defined 
with respect to the bare functions in the following way 



^ R(n) iPi(n),Qi{n), gn, fJ^n) — Z ^^^^Z ^2(^n)^ ' (Pi(n) , <5i{n) , Aon, /iOn, A„) . (53) 

We shall make reference to the cutoffs implicitly from now on. Although the regularization 
method utilizing cutoffs explained in the previous section is useful in many instances, we shall 
use dimensional regularization. Away from the critical dimension, the ultraviolet divergences 
of the original bare vertex parts can be represented by inverse powers (poles) in the variable 

in dimensional regularization. This manner of manipulating infinities will be shown in a 
moment to be particularly simple in the calculation of critical exponents. 
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The inductive proof of renormalizability of Eq.(53) at all orders in perturbation theory 
at the critical dimension (where the theory is renormalizable) was already demonstrated in 
the last section. Below the critical dimension, ai d — dc — e^, the theory is less divergent 
(superrenormalizablc) than its formulation at the critical dimension. Thus, we do not have 
to worry about the details of the rigorous proof of multiplicative renormalizability of the 
theory at o? = o?c — in analogy with our discussion in the previous section at d — dc, but 
take on this property as valid in that case as well. Instead, we shall give indirect evidence 
of this renormalizability by proving explicitly that the /3-functions and Wilson functions are 
all finite at the specific perturbative order we are interested in the present work. 

It is interesting to express the renormalized and bare coupling constants in terms of the 
dimensionless couplings Un, that is, = Un{nl^)^ ■> Aon = uonil^''^)^ ■ Moreover, we 
can write all the renormalization functions in terms of Un. 

The dimensionless bare coupling constants and the renormalization functions ^^(n), 
-^</>2(n) — 2^(i>{n)Z^'2(n) Can be expanded in terms of the dimensionless parameter Un up to 
two-loop order as 



These ingredients can be used to discuss the Callan-Symanzik-Lifshitz equations analogously 
to the analysis performed at the critical dimension. They shall provide a better comprehen- 
sion of the scaling behavior of the solutions in the anisotropic and isotropic sectors at the 
ultraviolet fixed points. 

A. Anisotropic 

There are at least two ways to derive the CSLE (at the and) away from the critical 
dimension d = 4 + X]n=2[ ''"n^'' ]^» ^^^^ manner, one can take a total derivative of 

the renormalized vertex part (including composite operators) with respect to the logarithm of 
the renormalized mass /x„ characterizing the m„-dimensional subspace at fixed bare coupling 
A„ and cutoff A„, in conjunction with the normalization conditions. The second trend is 
to take a derivative of an arbitrary bare vertex part with relation to the bare parameter 



(54a) 
(54b) 
(54c) 
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yUon) expressing everything solely in terms of renormalized quantities via normalization 
conditions, with bare coupling A„ and cutoff A„ kept fixed. The latter is almost identical 
to our procedure at the critical dimension dc = 4: + ^^^^[ ^"""^'' jm^ explained above. Both 
strategies lead to the corresponding Callan-Symanzik-Lifshitz equations: 

(A^n^ +/3n(fl'n,/in)^ - y7</.{n) + L7<^2(„) )rJJ[;^^ , Qi(„) , = (55) 

/"n"(2^ - 70(n))r}j(„) {Pi{n), Qt(n)\ 

These equations resemble the CSLE at the critical dimension, but now the renormalized 
coupling constants Qn are dimensionful. This provokes a discrepancy between the beta 
functions defined at the critical dimension and its corresponding version at 7^ 0. When 
expressed in terms of the renormalized parameters, the coefficients appearing in the CSLE 
equations satisfy the following relations: 

Pn {gn ,fJ'n) = /^n , (56) 

dlnZ^(n) f.-. 

l<t>{n)[9n, fJ-n) = fJ-n Tj , (57) 

7<^2(„)(5'„,/i„j = -fin — ^ • (58) 

Replacing the definition gn = ii„(/i^")^ into the above equation for Pnidn, f^n) we find 

d Ou d 

Pn{gn,IJ'n)^ = {fJ'n^)\0n^ h rieL5(„— , (59) 

dgn Ofln OUn Ogn 

where the derivative in the first term is calculated at fixed bare coupling constant Aon- It is 
equivalent to the dimensionless function /?„(«„) = (/inf^)Ao„- Now define the Gell-Mann- 
Low function for 7^ in terms of the dimensionless coupling according to 40, 4l| 

[Pn{gn, /in)]GL = -neign + /?n(fi'n., ^n) , (60) 

where the last term can be identified with the value of the function at the critical dimension. 
The solution of the Callan-Symanzik-Lifshitz equations for the vertex parts away from the 
critical dimension will possess a scaling limit provided it can be expressed entirely in terms 
of dimensionless quantities. We can get rid of the undesirable contribution coming from 
dimensionful couplings if we start from scratch with the Gell-Mann-Low function in Eq.(64). 
In that case, Eq.(59) turns out to be 

d d 

[(3n]GL{gn)-Q^ = PuM^;^. (61) 
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Consequently, the Callan-Symanzik-Lifshitz equations away from the critical dimension 
can be rewritten as 

(l^n-gj^ + PniUn)^ - y7</.(n) + L-f^'^{n))^^R(n){Pi{n), Qi{n),Un, Hn) = (62) 

Let us take a closer look at the dimensionless functions Pniun)- Since the derivatives 

of the rcnormalizcd masses in terms of dimensionless renormalized coupling constants are 
taken at fixed bare coupling constant, we can use the identity 

in order to trade partial derivatives. Using last equation, we find the following result for the 
fiow functions in terms of dimensionless parameters 

Pn{un) = -neL{^ ) . (64) 

Employing this result, the remaining renormalization functions can be rewritten in terms of 
dimensionless quantities in the form 

l,p{n){Un) = Pn (65a) 

l,p^n){Un) = -Pn ^^^^ (65b) 

l<f>\n)\Un) = -Pn • (65c) 

The equations (62) are the Callan-Symanzik-Lifshitz equations for anisotropic vertex 
parts describing generic competing systems of Lifshitz type with arbitrary composite op- 
erators which are multiplicatively renormalizable. Note that T^^^^ is precluded from this 
discussion. The independence of renormalization fiows in each subspace characterized by 
their corresponding bare(/renormalized) masses and coupling constants are manifest in this 
construction of anisotropic criticalities. 

The main obstacle to encounter a general solution of the CSL Eq.(62) is the appearance of 
the inhomogeneous term T^^^^~^^\pi(^n),Qi{n),0, dn, l^n) in its right hand side (rhs). It corre- 
sponds to the insertion at zero momentum in the vertex Tij^^^'' (Pi(n) > Qi{n),gn, l^n) due to the 
action of the derivative with respect to the mass parameter characterizing the n-th competing 
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subspace in the anisotropic cases. For the sake of simphcity, take L = 0. Generally speak- 
ing, this derivative produces an additional propagator in r^^^"*"^^ , (5i(n)5 0, yU„) when 
we compare it with T^^^^\pi(^n), Qi{n), 9n, l^n) ■ The limit oo describes the ultraviolet 

behavior of such vertex parts. Thus, r^;^^^(pi(„), 0, gf^, /i„) is of order Pj^^"r^^^(pi( 
for large p„, up to powers of Inpn, at every order in perturbation theory. We make the 
assumption that these logarithms do not add up to compensate the factor p^^^". 

All the momenta involved in a given Feynman diagram are Euclidean. Taking a fixed 
nonvanishing Euclidean momentum ki(^n) , we can reach the ultraviolet region of any diagram 
by performing a scale transformation like = Pnki(n), where p„ is a dimensionless flow 
variable in the limit p„ ^ oo. Interesting anisotropic configurations have the following 
structure: there are n independent subsets of momenta which go to infinity independently. 
Since the n independent mass subspaces have a well defined zero mass limit and all the 
momenta appearing in r^^-|(pj(„), are nonexceptional in the ultraviolet regime (except 

for the zero momentum of the inserted ^operator in the vertex part {Pi{n) , 0, gn, f^n)), 
we can apply the Weinberg's theorem 42] in each subspace, which guarantees that the rhs 
of Eq.(62) can be neglected order by order in the perturbative expansion. 

The limit p„ cxd implies to take all internal momenta in a given diagram as large as 
possible. This goal can be scored by regulating the integral with cutoffs A„ in the limit 
A„ —>■ oo. Hence, the integral becomes a homogeneous function of the mass The regions 
in momentum space where the rhs in Eq. (62) can be neglected when compared to the 
left-hand side (Ihs) come from the limit oo which is identical to the ultraviolet limit. 

This is indeed the case, as we have already shown in Sec. Ill using cutoffs. Moreover, it 
will be demonstrated in the Appendixes via dimensional regularization that the ultraviolet 
divergences are associated to this region in momenta space where the loop integrals shall be 
computed. 

Recall that several dimensional redefinitions of the momentum components along the 
various competing subspaces have been performed as discussed in Sec. II (see also |36|). If d 
is the spatial dimension of the system, these redefinitions along with the Lifshitz conditions 
translate themselves into an "effective" space dimension for the anisotropic situations of 
generic competing systems, i.e., d^fj = d — X]n=2[^^]"^n- Elementary dimensional analysis 
implies that, under scaling, the vertex parts possess the following behavior 
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il(r)vKnft'i(n)) "n; Pnj — Pn '- R{n)\'^i{n) i '^n\Pn ) i )■ 

Pn 

The solution of the asymptotic part associated to the vertex function satisfying the ho- 
mogeneous CSL equations has the property 

(AT) u N du' 

KsR(n)^Hn),Un,^) = exp[-— / 1 4>{n){u'M) ITTZF^^ (66) 
Pn ^ Jun h'n\"'n) 

il(n)(^«(")' '^n{Pn) , Pn) , 

where 

Prom now on we omit the subscript for asymptotic in the vertex parts satisfying the homo- 
geneous Callan-Symanzik-Lifshitz equations. The eigenvalue conditions Pn{unoo) = yields 

the ultraviolet nontrivial fixed point. Exactly at this value of the coupling, which is indepen- 
dent of the competing subspace under consideration in the anisotropic cases, the solutions 
of the CSL Eq.(62) under a scale in the external momenta can be written as 

(n h-, -.11 II ) - ^"^^^+^'* 2.„=2i— 2—" ) 2 (CO) 

"'^il(n) (^*(")' ^nco-i Pn)- 

The dimension of the field can be defined by 

^'^{Pnhin)) = p::[('^-^-[^]-)-^'^^(")]r|::;(fc.). (69) 

Due to the interactions, the field develop an anomalous term, which can also be defined 
through the relation d^^^i) = '^~^"=^j'^''"" _ i + 2^. Comparing Eqs. (68) and (69), one 
can easily verify that the anomalous dimension r]n can be identified with 70(n)(iinoo)) i-©-) 

^7" 0'</i(n) (^noo) • 

Plus, the anomalous dimension of the composite operator 0^ can be extracted by analyz- 
ing the vertex parts including composite operators. The asymptotic behavior of the vertex 
parts at the fixed point w„ = Unoo coming from the solution of the CSL Eq.(62) has the 
following simple scahng property {{N, L) ^ (0, 2)) 



R{n) yPni^t^nji PnPi(n)i '^nooi PnJ — Pn 

X Pn ji{n)\i^i{n)-iPi{n)-i'^noo-i Pn)- \'^) 
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If we write the coefficient in the rhs as pn "=2 n " « ^ (n) ^ discover that 
d(f)2(n) = — 2n + 7(^2(„)(m„oo). The correlation length exponents can be computed through the 
identification z/"^ = -(i^2(„) = 2n - 7<^2(n)(wnoo)- 

This method with several independent mass scales generalizes, in a nontrivial way, the 
previous method developed for m-axial anisotropic Lifshitz points where only two mass 
scales are present. 



B. Isotropic 

The existence of only one type of competing axis in the isotropic criticalities make the 
analysis simpler for the cases d = rUn- As before, we shall briefiy discuss this case by focusing 
on the analogy with the competing sector labeled by n of the anisotropic treatment already 
discussed above. We start with the mass scale dimensionless coupling constant Un and 
subscript n in all vertex functions. The critical dimension is An, whereas the expansion 
parameter is eL = An — d. The effective space dimension is At the (eigenvalue condition 
f3n{unoo) = 0/) fixed poiut UnoQi a scale transformation in the external momenta implies that 
the vertex functions present the following property: 

The dimension of the field is defined by 

r|:,'(p„fe.,„,)=p:«"'-™*«'rW(*.). (72) 

Then, it follows that = ^ — 1 + ^, which in turn implies that the anomalous dimension 
of the field (rjn = i]Ln) satisfies the identity rjn = 7</,(n)(w„oo)- 

Now, let us consider the vertex parts including composite fields. The identification of 
the anomalous dimension of the composite operator 0^ can be done through the following 
steps: at the fixed point w„ = Unca the asymptotic behavior of the vertex parts under scale 
transformation can be written as 

r(^'^)fo a) - "(^+^-^^)-^^)- "^^^n'""°°V ^7,2,„,K^) 

R{n) \Pn'^i{n)'> HnPi{n)'> "'nooj hl'n) — Hn K'^) 
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Then, write the coefficient in the rhs as pj " 4>^(n) ^ learns that (i02(„) = 

—2n + 7(/,2(„)(-u„oo) and the correlation length exponent (z/„ = u^n) in the isotropic case can 
be identified with (i<^2(„) through u'^ = -(i<^2(„) = 2n - 702(„)(m„oo)- 

Finally, the beta function does not have the global factor of n as in the competi- 
tion directions of the anisotropic case. Rather, in terms of dimensionless parameters 

These are all informations required for the calculations of the critical exponents by dia- 
grammatic means. This theme will be developed in the following remaining sections. 

VI. CRITICAL EXPONENTS IN THE ANISOTROPIC CASES 

In order to calculate the critical exponents for anisotropic critical behaviors of arbitrary 
competing systems, we shall use the results in Appendix A for the massive Feynman graphs 



using the generalized orthogonal approximation. As already demonstrated in 35|, l36| , this 
approximation is the most general one consistent with homogeneity in the several indepen- 
dent scales of external momenta. In the massive case introduced in the present work, the 
orthogonal approximation is also consistent with homogeneity in the various mass scales 
associated to each competing subspace. We shall employ the normalization conditions setup 
of Sec. II in our computational procedure. 

To begin with, we use Eqs.(3) and (53) to fixing the normalization functions 2'^(n), Z^2(^n) 
in powers of u„ in conjumination with the appropriate Feynman diagrams at the loop order 
needed. In the expressions of the one- two and three-loop level of the integrals I2, 1'^, li^ 
given in Appendix A, we emphasize that a geometric angular factor appears every time we 
perform a loop integral. It is explicitly given by [5(d_j^^^^„„)r(2-^„^2 ^)0^n=2 )] 
and can be absorbed in a redefinition of the coupling constants. We can perform those 



redefinitions in order to discard it from our considerations henceforth. We then obtain the 
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renormalization functions in terms of those loop integrals in the following form 

+ .[(AT + 8)72]' 

^ (7V2 + 67V + 20)/| , (5iV + 22)/4 , ^j,. 
7 1 , {N + 2)I' (7V + 2)(7V + 8)(/24-|) 

^<^(n) = 1 + Un + U^, (74b) 

- (iV + 2)/2 „, 

^</)2(n) = -I- H g Un 

^^ (A.^ + 77V ^ 10)7. _g+2)^(j^^^^^j^. p^^, 

In dimensional regularization, the coefficients of the various terms in powers of Un present 
poles in e^. They cancel in the evaluation of /3„ and the critical exponents. When we 
combine Eqs.(54), (64) and (65) together, we obtain the and Wilson functions in every 
subspace. We can write them through the following expressions 

I3n = -neLUn[l - ainUn + '^{al^ - a2nWn], (75a) 
7</.(n) = -neLUn[2h2nUn + (S^sn " 2h2nain)ul\, (75b) 

Icpin) = neLUn[cin + (2C2n " Ci„ - ainCln)Un]- (75c) 

Using the explicit values of the integrals presented in Appendix A, we can determine the 
above mentioned coefficients by using the normalization conditions as functions of those 
integrals calculated at zero external momenta. When the results in Appendix A are combined 
with Eqs.(54) and (74), we conclude therefore that 

am = ^^^^[1 + {hm, - l)eL], (76a) 

_ (jV + 2)(7V + 8) (jV + 2)(jV + 8) 1 13 
129fe| + lO&I ' 4''"" + 48^' 

Cm = ^[1 + (/im, - 1)6^], (76e) 

(Ar + 2)(Ar + 5) (Ar + 2)(Ar + 5) (4Ar^ + 25Ar + 34) 

- 364 18^^ ^ "^^ " ^ 72^^ ■ ^ ^ 

Just as in the case of noncompeting and m-axial universality classes, the integral I encoun- 
tered in Appendix Eq.(A9) appears explicitly in the coefficient 62n Eq. (76c). We emphasize 
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that these expressions differ from Eqs.(67) of Ref. [36l| obtained in the massless theory at 
nonvanishing external momenta in each subspace. Since the normahzation conditions in 
the present work are defined at zero external momenta as well as nonvanishing mass in the 
specific subspace labeled by n and recalling that these coefficients depend upon the nor- 
malization conditions, this disagreement in their values using either renormalization scheme 
should not be surprising. 

Inserting the values of the coefficients given in Eqs. (76) into (75a), /?„, can be written as 







_P^„JH0(.4). (77) 

Using last equation, the eigenvalue conditions Pniunoo) = lead to two solutions: a trivial 
zero as well as a nontrivial zero of order of each j3n characterizing the independent 
noncompeting (n = 1) and competing (n = 2, L) subspaces. Remarkably, they correspond 
to the same value of the coupling constant (uioo = ^200 = ••• = mloo = Unoo), i-e.. 



1N (9A^ + 42) 
-{h„,^ - 1) + ' 



(78) 



(8 + Ny 

Substitution of this value in the functions 7<^(„) and 7^2 („) together with the coefficients 
given in Eqs. (76) allows us to find the critical exponents rjn and Un- Indeed proceeding as 
indicated, we first find 

, Ti 2 N + 2 . ,6(3A^+ 14) 1,^ , , 



Note that these expressions are exactly the same of those obtained previously in [36|, Eq.(73) 
therein, upon the identification 70(„)(u„oo) = Vn- 

The exponents can be encountered from the expression of the anomalous dimension 
of the composite operator 0^. As we know, 7(^2(„) = 7(^2(„) + 7<^(„). Thus, it follows that 
-d(i>^{n) = 2n - 7<^2(„)(u„oo) - l^{n){unoo)- Usiug Eqs.(76) once more, the following 



-1 



intermediate result can be verified: 

7<^2(n)(Mn) = n ^^'^'^\ n[l + eiihrriL " 1) " (80) 

Using the value Unoo into the last expression along with 70(n)(M„oo) from Eq.(79), we get 

_ 1 , iN + 2) , 1 (iV + 2)(iV^ + 23iV + 60) , 
" 2n ^ 4n(iV + 8) ^ 8n (iV + 8)3 ^""'^ 
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These exponents are in exact agreement with those calculated previously using a massless 



theory 36|]. Note that the properties of strong anisotropic scaling 43| rjn = nrji and z/„ = 
(^)z/i already encountered in the massless case are reproduced using quite a different method 
and generalizes the more restricted m-axial Lifshitz universality class 27|, ImI] • 



VII. ISOTROPIC CRITICAL EXPONENTS IN THE ORTHOGONAL APPROX- 
IMATION 

Although the isotropic diagrams can be calculated exactly without the need of any sort 
of approximation, in this section we shall evaluate the isotropic exponents for two reasons. 
Firstly, for completeness: pursuing the analogy with the anisotropic cases within an ap- 
proximation which keeps the desirable concept of homogeneity in the mass(es). Secondly, 
simplicity is a key ingredient, since possible applications in other field theories with higher 
derivatives might be worthwhile using this simple setting of computation. Due to nontriv- 
ial features of the exact computation, we postpone the discussion of the calculation of the 
integrals and of the exact critical exponents to the next sections. 

We now have only one type of subspace to be integrated over, along with only one type 
of mass parameter for each isotropic behavior labeled by the number of neighbors n coupled 



via competing interactions according to the CECI model introduced in [35|] to understanding 
the most general magnetic systems with competing interactions. We can take advantage of 
the discussion in the previous section: we focus on the nth subspace taking into account 
of course the particularities of the isotropic renormalization functions as pointed out in 
Sec. V. In addition, we shall use the results in Appendix B for the Feynman integrals in 
order to compute the exponents rjn and z/„ via perturbative expansion. Each loop integral 
originates an angular factor of Sm„, the area of an m„-dimensional unit sphere, which shall 
be absorbed in a redefinition of the coupling constant. In order to make a clear distinction 
among the various isotropic behaviors and the nth competing subspace of the anisotropic 
behaviors, we shall replace the perturbative parameter in the isotropic case, i.e., e„. 
Furthermore, this slight change of notation shall be useful to perform a comparison with the 



critical exponents already evaluated in the massless theory in Ref. 361]. 



First we list the analogues of Eqs.(75) for the renormalization functions in the isotropic 
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cases, namely 



4], 



lcf>{n) = -e„M„[262„M„ + (363„ - 2b2nain)ul], 
7</.2(„) = e„M„[cin + (2c2n - Ci„ - ai„Cin)M„]. 



(82a) 
(82b) 
(82c) 



Employing the results of the previous section along with the output obtained in Appendix B 
utilizing the orthogonal approximation in the computation of the appropriate loop integrals 
for the isotropic behaviors, the various coefficients can be expressed in the form 



Cln 
C2n 



A^ + 8 
,A^ + 8, 



6e„ 

(AT + 2) 



,2iV2 ^41A^ + 170, 



144ne.^ 
(Ar + 2)(Ar + 



+ 



_ (iV + 2)^ 
144n2 ' 

7(A^ + 2)(iV + 
5184n2e„ 



- i.„]. 

6e„ ^ 2n ^' 
(A^ + 2)(A^ + 5) (A^ + 2)(2iV + 13) 



(83a) 
(83b) 
(83c) 
(83d) 
(83e) 
(83f) 



36e^ 72ne„ 

Note that these expressions are different from Eqs.(90) derived in Ref. 36|] using the orthog- 
onal approximation in the massless limit. Since e„ = 4n — and substituting the above 
results in the expression for we find 



/3n 



(AT^ 

(3iV+W) ' 
+ 12n 



■(1 - 7r^n)Un 

In 
OK). 



(84) 



The eigenvalue condition (3n{unoo) = permits us to find out the nontrivial fixed point of 
the dimensionless coupling constant, whose value is given by 



- 6 I e„ 



1 {m + 42) 

2^ (8 + Ar)2 



(85) 



Replacing this fixed point in the expression for 70(n) together with Eqs.(83), we then have 



7<^(n)l 



Mr: 



2^V + 2_ 1^ 6(3Ar+14) 1^^ 



2n{N + 8)2 ' 



■n' (Ar + 8)2 4^ 



(86) 
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This value corresponds exactly to the exponent rin previously obtained using the renormal- 



ization group equation in the massless theory 
(83) in the definition of 7(^2(„), we obtain 



35, 



36|. Furthermore, using again the results 



Replacing Eq.(85) into last equation, we get to the following result 

N iN + 2) 6(A^ + 3) , 

%Hn)M = Jf^fnil + ^^^^^M (88) 

Using the relation i/„ = {2n — 'y^2(^n){unoo) —1<i>{n){unoo))^^, we find out the correlation length 
critical exponent 

1 , iN + 2) , 1 (iV + 2)(iV^ + 23iV + 60) ^, 

The expression for 702(„)(m„oo) in Eq.(88) is the same as the one associated to a scalar theory 
in the massless limit, computed at the fixed point using the renormalization group equation. 
Besides, Eq.(89) corresponds to Eg. (95) for this exponent in the orthogonal approximation 



using the massless method of Ref. [36 1. 



We can compare the results given in Eqs.(86) and (89) of the exponents rjn and Un 
for generic n with the previous massive method recently introduced to treat the n = 2 
isotropic case corresponding to m-axial Lifshitz points in Ref. [3ll] employing the orthogonal 
approximation. Two misprints in the expressions for the critical exponents 772 = Vla, and 
1/2 = namely Eqs.(95) and (98) in [31|], respectively, took place in that paper. A wrong 
factor of 2 occurred in Eq.(95) as the coefficient of the dependent fraction of the 0(e|) 
contribution. In Eq.(98) an incorrect factor of | also appeared in the 0(e|) contribution 
therein. We emphasize that the correct expressions for those exponents are given by Eqs.(86) 



and (89) for n = 2 with in [311] identified with €2 herein. 



VIII. ISOTROPIC INTEGRALS IN THE EXACT CALCULATION 

In the present section we shall discuss and compute the isotropic integrals for arbitrary n 
in order to calculate the critical exponents in Sec. IX. This computation in the present massive 
method is by far much more complicated than its counterpart in the massless framework of 
Ref. [36| . Since this setting for the case n = 2 was not explicitly demonstrated in our previous 
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work 



3l[ | (the solution to the integrals were only quoted in one of the appendixes to that 



paper), we take this opportunity to fill this gap by working out more complicated cases 
along the same lines of reasoning. Here we shall describe in great detail all the technicalities 
involved in the evaluation of these integrals. 

As we are going to see, the two- and three-loop contributions for the two-point functions 
represented by the integrals l'^ and J5, respectively, present no problem in their calculation. 
In fact, they can be solved in terms of the expected poles and regular terms accompanied by 
multiparametric integrals which do not need to be computed explicitly, but cancell out in the 
renormalization algorithm. The issue is the computation of one- and two-loop contributions 
to the (four-point vertex part) coupling constant. 

We start by considering the integral J2. Since it always appears as a subdiagram in 
higher loop contributions of arbitrary IP/ vertex parts, we first attempt to compute it at 
nonvanishing external momenta. In this example we shall have a precise idea of the difficulty 
in this computation. 

At nonvanishing external momenta K, I2 is given by 

((A; + ir)2)« + l)p2)- + l] ' ^ ^ 

From elementary complex algebra, we can use the identity (/c^)"-|-l = {k'^—ri){k'^—r2)...{k'^ — 
r„), where r/ is the Ith complex root of the equation (A;^)" -|- 1 = 0. 

To see this procedure at work, we apply it to the simplest nontrivial case which occurs 
for n = 2. In that case, using Feynman parameters, the propagator can be written as 

r(2) TIT^, (91) 



(A;2 — ri)(A;2 — r2) Jq {k'^ + m'^y 

where ml = (r2 — ri)x — r2. In particular, using Feynman parameters to fold denominators 
with different roots amounts to write the last expression for ml in different ways as functions 
of Vi and r2. Analogous remarks are valid for multiple roots. The final answer for the specific 
integral we are interested in is naturally independent of these maneuvers. For the case n = 3, 
a similar reasoning leads to the following representation for the propagator 

r(3) [ [ dxidx2X2 ^ .0 , (92) 



(/c2 _ ri)(/c2 _ r2)(A;2 _ rg) Jq Jq {k^ + ml, 

with m^^ = {m^^ + t^)x2 — T3 and so on. For arbitrary positive integer n, we find by the 
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same token the following result: 

_ ^ = T{n) [ ... [ dxidx2X2...dxn-ixlzljYY-——2 v^' (93) 

where = {ml^_^+rn)xn-i—rn- Now let us insert Eq.(93) inside Eq.(90). Consequently 

I2 becomes 



I2 = T{n) / ... / dxidx2X2...dxn-ixl_^dyidy2y2---dyn-iyn~i 
Jo Jo 

X 1 d"^"k ^ (^A) 

We then utilize another Feynman parameter, say t, in order to fold the two denominators. 
Integrating over the loop momenta 



36| . we get an angular factor Sm„ which can be dis- 
carded/omitted through the redefinition of the coupling constant just as before. The reader 
should be warned that we shall use this fact in all loop integrals henceforward. Performing 
the e„-expansion of the prefactors (which are simple Gamma functions), I2 turns out to be 

h= [I + - %l}{2n))\ j ... j dxidx2X2...dXn-ixlz\dyidy2y2---dyn-iyn-i 

^ Jo Jo 

X dt[t{l - t)r'[t{l - t)K' + {ml_^ - mlJt + mlJ^, (95) 

where ipiz) = '^^"j^^^-' is the digamma function. One could think that matters get simple to 
grasp if we set if = (which is actually the symmetry point we shall be concerned in this 
massive setting), namely 

-{ipil) - ilj{2n))] ... I dxidx2X2...dXn-ixlzldyidy2y2---dy 



r(2n) 



fcn ^ JO JO 

X dt[t{l - t)r\{ml_^ - ml^Jt + m^„_J^, (96) 

but this is not so. The reason for the difficulty even for = is that the complex roots 
of the polynomial in k corresponding to the propagator depends explicitly on the value of 
n: we can not just go on, within this method, without specifying the value of n. This is the 
main obstruction to find the result of the integral in terms of a pole in together with a 
simple regular term for general n. 

The simplest way to find the values of the integral for arbitrary n is trying to discover 
a recurrence formula for this integral by analyzing the cases with fixed n at the above 
symmetry point. Two cases are already known. The case n = 1 which corresponds to the 
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standard quadratic 0^ theory, whose integral is given by I2SP = ^(1 ~ The case n = 2, 
previously derived in Ref.|31i is given by I2SP = ~(1 ~ It would be interesting to 
compute this integral for higher values of n in order to see if the method can be trusted for 
arbitrary n, provided the value of n is fixed. 

With this idea in mind, let us compute the integral for n = 3 in the first place. The 
three complex roots are ri = —1, r2 = ^ + = A and = A*. We find useful to melt 
firstly the contributions of r2, r-^ and folding the resulting expression with ri term afterward. 
Expansion of the digamma functions in the prefactors, we find 



120 137 

hsp = [1 ~ TT^TT^s] / '^^ '^^ / dww I dt\t{l—t)] 



,-•3 11 ii.'.r, I ii.z. I (1,11 II I (1,111 III I (;,/,!/, I I — /, 1^^ 

es ' 120 

X [{^V3x - ^ - ^^)y{l - t) + {^V3z - ^ - + 1]^ . (97) 

When we perform the parametric integrals, we set €3 = in the powers of the resulting 
complex numbers which remain. The above expression is thus equal to 

hsp = -(1 - (98) 

63 

For good measure, let us analyze the case n = 4. Since this case is more involved and the 
calculations are rather long (though straightforward), we begin with the most basic facts. 
The complex roots are ri = | + ^-^ = B, r2 = —B*, = —B and r4 = B*. The propagator 
can be written as 

where ml = {[{2x — 1)B* — B]y + 2B)z — B. After using the Feynman parameters, the 
integral I2SP then becomes 

/"I r dJ^'^k 

TiAY I ^^^ydrrd^^'dss' j ^ ^ (100) 

with ml = {[{2w — 1)B* — B]r-\- 2B)s — B. Using another Feynman parameter t, integrating 
over the momenta and absorbing the angular factor 5*^4 in the redefinition of the coupling 
constant, we obtain 

rf8) 363 

hsp = [1 £4] / ••• / dxdwdyydrrdzz^dss'^dt\t{l — tyf 

£4 280 Jo Jo 

X [{ml - ml)t + m^^^ . (101) 

We are thus left with the task of calculating a large number of elementary integrals, which 
makes the whole process a tedious one. Obviously, the order of the integrations chosen was 
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'25P 



the following: integrate over w, y, r, z, s and t. In the end of the day, set 64 = in the 
powers of the remaining complex numbers to get to the result 

/25P=^(l-§). (102) 
£4 o 

From the cases n = 1,2,3,4, we discover that the method is reliable for arbitrary high 
values of n. There is no change of pattern in the calculations for different higher values of n. 
Perhaps the only undesirable feature is the proliferation of elementary parametric integrals 
whose number increases with increasing n. Therefore, for arbitrary n we conclude that 

W=i(l-|i). (103) 

We then find out that the orthogonal approximation discussed in Appendix B is exact at 
one-loop level for the massive method. 

Let us turn our attention to the integral I3. In a preliminary stage, we leave the external 
momenta arbitrary and in that case it is given by: 



(104) 



X dz^dz2Z2...dz„.lz:zldt[t{l - i)]"-^-^ / , , ,,,, , (105) 



iiih + k2 + KYY + 1) [{klY + l\[{klY + 1] 
Now, performing the integral over the loop momenta k2 by conjugating the results expressed 
in Eqs.(93) and (95), it is easy to see that 

h = [1 4_ ^(^(1) - ^(2n))] f ... f dxidx2X2...dxn-ixlzldyidy2y2---dyn-iynZl 

^ Jo Jo 

d'^^ki 
where 

"f^f — 7-. — ^ 

t{i-t) 

and ni\^_^ = + Tn)zn-i — Tn- Wc then use another Feynman parameter, say u, in 

order to integrate over the loop momenta ki to find 
r(2n)r(2n-f)r(-n + e„) 

en 
2 

X dyxdy2y2---dyn-iyZ^ldzidz2Z2...dzn-izlZldt[t{l - t)X'~'^~'^ 
X duu'^-^{l-uf-\u{l-u)K'^ + {ml-ml^_^)u + ml^_^f-'-^ . (106) 



r(2n)r(2n- ^)r(-n + e„), ,,,,,, /"^ 
^3 = orr^n^ ^[1 + f (^(l)-^(2n))] / dx,dx2X2 dx^-.x^l 



From this expression we can compute the derivative at zero external momenta, i.e., /g = 
QK^n \K'^ri=Q. First, we expand the F-functions in the prefactor using the appropriate identities 
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for these functions. Second, after taking the derivative and setting the external momenta 
to zero, the first term in the last bracket contributes to the integral over u, whereas the 
powers of the remaining parameters are expanded in e„ through the elementary identity 



a ^"{u,t, Zn-i) = 1 — €n In a{u,t, Zn^i). Thirdly, the identity |44| 



/ dxx^'-\l-x'Y~Hnx = -Bi^,v)Wji^)-^lji^ + v), (107) 
Jo j,l J, ^. J. 

where B(x,y) is the Euler beta function, will be useful to our purposes. It allows to separate 
the pole and a regular term by solving all the parametric integrals together with a regular 
term which depends on only one multiparametric integral. Collecting together this set of 
steps, J3 turns out to be 

(-irr(2,,)^ -3 'fix '^l 

p—o p— 1 

where the remaining multiparametric integral H' is given by 

H' = dxidx2X2 dXn-ix'^Zldyidy2y2---dyn-iynZldzidz2Z2...dZn-izllzldt[t{l -t)]""'^ 

Jo 

X duu-'\l - ur-Hn{[^^:<-^^^ . (109) 



The above integral is the generalization of the integral H Eq.(C5) from Ref. 31i appearing 
in the exact computation of the analogous two-point diagram for m-axial Lifshitz points. 

After this discussion, it is straightforward to perform the three-loop integral contributing 
to the two-point function, namely 

d"'"kid"'"k2d"'"k3 

(((fci + k2 + KfY + 1) (((A;i + fcs + KfY + 1)((A;?)" + l){{klY + l)((A;i)" + 1) 

(110) 

Note that the internal subdiagram of the four-point function appears quadratically in the 
above integral. Using Eq.(93) in conjunction with Eq.(95), we get to the following interme- 
diate result 

T{2n)T{n) 



L. 
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-[1 + e„(?/'(l) - ?/'(2n))] / ... / dxidx2X2...dxn-ix'^_\dyidy2y2---dy 

Jo Jo 
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Proceeding as before, integrate over the remaining loop momenta utilizing another Feynman 
parameter. This leads to 

r(2n)r(2n - f )T{-n + ^) . 



[l + e„(V'(l) - V'(2n))] / dxidx2X2 dxn-ixl_\ 

Jo 



>< dyidy2y2---dyn-iyn-ldzidz2Z2...dzn-izl_ldt[t{l-t)f ^ 



X duu'^-\l-uf-\u{l-u)K^ + {ml-ml^_;)u + ml^_;\''-^-^ . (112) 

Wc then compute /g = gf^lA'an^o- Taking the derivative, performing the expansion of the 
F-functions, using it^" = 1 + Cnlnu along with a similar expansion for the brackets (naturally 
taken at = 0) and employing the identity (107), the final result for the integral is 

3n— 1 T Ji— 1 2n— 1 



-l)"F(2n) 
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3F(n+l)F(3n)e2[^ + ^"( ^ I] p + 2^ 2p + 5Z 2p^^ 

^ ' ^ ' " p=3 p=l p=l 

31 (n + l)e„ 

where H' is given by eq.(109). 

Finally, let us compute 1^ at zero external momenta, namely 

d^-hd^-k2 

[{klY + iY{{klY + i)[{{h + k2fY + i) ■ ^ ^ 

The same kind of limitation taking place for I2 occurs for I^: the lack of a solution for general 
n due to the fact that the roots of the polynomial in the momenta in our new representation 
of the propagator depends explicitly on n. To see this, integrate first over the loop momenta 
k2 just as explained above, that is, use Feynman parameters to rewrite 1^ as 

h = ^^^^^^^[1 + ^(V'(1)-V'(2n))] /'... f dx,dx2X2...dxn-,xlz\dy^dy2y2...dyn-iyT-l 

^ Jo Jo 

.Mlt(l-tT-'-*l^^^^^,^^^^^,^. (115) 

If we use new sets of Feynman parameters Zi to work out the first propagator, this operation 
produces the expression 

' - l^(^''^l^H[i + |(^(i) _ ^(2n))] 2' ... dx^dx2X2...dXn-lx:zldy^dy2y2...dyn-ly:,Zl 



X dt[t{l - i)]"-^-^dZi^i(l - Zi)dZ2zl{l - Z2)...dZn-lzl'rf{l - Z^-l) 

f d'^^kx 



(116) 
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Performing another integral over the remaining loop momenta the integral can be entirely 
represented by parametric integrals, that is 

^ r(2n)r(2n-f)r(6„) e , m ^^^ f\ , , 

[1 + —{ipil) - i'{2n))\ / dxidx2X2 



2r(f)e„ ' 2 ^0 

X dyidy2y2---dyn-iynZldZiZi{l - Zi)dZ2zl{l - Z2).:dZn-lz'^Si{l - Zn-l) 

X dt[t{l - t)r'~'^duu'f-\l - ur-\{ml - mlJu + m^^^J— . (117) 

A standard procedure is to compute the integrand at m = 0: by summing and subtracting 
the integrand at this value, the difference between the integral at n 7^ and at m = is 
higher order in e„ and can be neglected 22, |45|]. We emphasize the factor Tn1^_^ depends 
upon the variables Zi. Therefore, we can integrate over the other parameters and performing 
the e„-expansion of the F-functions, we obtain the following expression: 

h = ^|^[l + 6„(^^(l)-i^(2n)-^(n))] 

dziZi{l-z{)dz2zl{l-Z2)..Azn-iz''^Si\l-Zn-{){ml^_X"'' ■ (118) 

In order to proceed from this point, we must specify again the values of n in order to 
determine the roots explicitly in order to solve the remaining parametric integrals. We 
already know this integral for n = 1,2, namely 

/.^^(l-f). (120) 

The cases n = 3, 4 can be worked out just as we did in the one-loop case. With the resources 
furnished so far, it is not difficult to obtain the following results 

/4 = A(1-^), (121) 
2ep 120^' ^ ^ 

/4 = ^(1-^). (122) 
2el^ 840 ^ ^ ^ 

The value of 1^ for general n which correctly reduces to the above particular cases is then 
given by 

h = ^{l + e^{D{n)-\)). (123) 

where D{n) = |^(1) - ^(n) + l^/'(2ra). 

This completes our task of calculating explicitly the exact isotropic integrals in the mas- 
sive case |46| necessary to compute the critical exponents. This aim shall be tackled in the 
next section. 
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IX. ISOTROPIC EXPONENTS IN THE EXACT COMPUTATION 



Now we apply the normalization conditions of the massive theory in pretty much the same 
way as worked out in the anisotropic and isotropic cases. In other words, the algorithm to 
determine the critical exponents is the same, but now we have to replace the values of the 
integrals calculated in last section. 

The coefficients of the several rcnormalization functions can be easily found. Indeed, 
using the Wilson functions appropriate to the isotropic cases the cofficients have the following 
expressions: 

,iV + 8,2 ,N^ + 26N + 108, 5N + 22A 

= '-fcT' " ' ' + - ^("') 

r(2n)'(^ + 2) 

^ ^ 36r(n+ l)r(3n)e„' ^ ^ 

= (-1)- r(2n)-(^ + 2) 3 ^ g Ij'f 1 _ r(3„)^,)|_ (12^) 



r(2n)2(iV + 2)(iV + 8) ^ ^ ^"-^ 

53n = (-1) 



^ n— 1 2n— 1 ^ 

648r(3n)62 + '"^"4 + |J ^ " |J p 

3n-l 
p=3 

In practice these results determine the several renormalization functions. They are given 
by: 

PniUn) = -uJen (l " —)Un + ( D(n) + -1 " ^ , \ , \ ) 

"L" 6 ^ 2n^ ^ 9 V ; V y I8r(n+ l)r(3n)^ 
X ul], (125a) 

+^(^ E ^ + E ^ - ^ nZn)H')u.] . (125b) 

p=2 p=3 ^ 

%^n){Un) = '^^^ ^n[l ^ ^ ^ -D(n)li„] . (I25c) 
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The fixed points Unoo are obtained from tlie eigenvalue conditions Pn{unoo) = 0, wliicli yield 



6 



It should be pointed out that neither the coefficients/renormalization functions nor the 
fixed points are equal to those obtained in Sec. VII of Ref.[36| for the massless theory: the 
normalization conditions are different and the difference is due to the nonuniversal feature 
of these functions away from the fixed points. Incidentally, we make the observation that 
even the fixed points are not universal but actually vary with the renormalization scheme 
employed. 

Replacing the fixed points back into 7(/,(ra)(Mn), the result is just the anomalous dimension 
Vn = 70(n)(Mnoo), which is glvcu by 

„ ( (iV + 2)r(2n)^ {N + 2)T {2ny F {N , n) 

^ ^ (Ar + 8)2r(n+ l)r(3n) " ^ ' (A^ + 8)2r(n + l)r(3n) " ' ^ ' 



where 



2n— 1 ^ n— 1 ^ 3n— 1 ^ 



p=i ^ p=i ^ p=i ^ 



Comparing this expression with Eq.(lll) from Ref. ,36|, we detected a misprint therein: 
an extra factor of ^ appearing there should be disconsidered. The correct expression for 
F{N ^ n) is given by the last equation shown above. 

Note that the integral H' cancelled out in the computation of rjn- Finally let us compute 
the exponent z/„. First, we write 7(/,2(„)(u„) at the fixed point, and then use the scaling 
relation among t'n, 7(/)2(„)('U„oo) and rjn to get the following result: 



2n 4n2(iV + 8) 4n2(Ar + 8)3 ' ' 'r(n + l)r(3n) 

^(iV + 2)(iV + 8)^ (129) 
2n 

3611 in the zero mass 



These universal quantities are the same as those obtained in Ref. 
approach. The other critical exponents can be obtained by scaling relations from the two 
exponents above and are also the same as those calculated from the zero mass setting. A 
rather interesting point is the following: while the exponents in the zero mass limit are 
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obtained analyzing the behavior of the vertex parts in the infrared, the exponents here are 
obtained in the ultraviolet regime since the fixed points m„oo are ultraviolet fixed points. We 
thus have complete equivalence of the zero mass theory renormalized at nonzero external 
momenta with the massive theory renormalized at vanishing external momenta. Universality 
is obeyed as expected. 



X. FINAL COMMENTS 



The introduction of a massive method to calculating critical exponents of generic com- 
peting systems is another step forward to a better comprehension of Lifshitz criticalities 
for several reasons. First, the massive method along with the Callan-Symanzik-Lifshitz 
equations is appropriate to proving the multiplicative renormalizability to all orders in per- 
turbation theory as demonstrated herein. This proof for generic competing systems is a 



nontrivial generalization of that already performed for m-axial Lifshitz points [3ll]- Besides, 
the massive framework with its complete equivalence to the zero mass treatment reflected 
in identical critical exponents with those in Refs. 36|] whether we use approximations for 
the appropriate cases or approach the calculations exactly is another exact manifestation 
of universality. Let us discuss some properties of anisotropic and isotropic cases separately 
and how they can show up in other field-theoretic contexts. 

It is worthy to emphasize some especial features that occur in both massless and massive 
renormalized perturbation theories constructed to describe the generic higher character Lif- 
shitz universality classes. As we have discussed, the computations in the massive theory is 
much more elaborate than in the massless method. Even though the Wilson functions de- 
termining the renormalized theory are explicitly different in both cases, they converge to the 
same value at the respective fixed points, which make them responsible for the same values 
of the critical exponents (and all universal quantities) in either scheme. This is not the only 
difference: while the massless IP/ vertex parts are scale invariant at the infrared attractive 
fixed point m* , their counterparts in the massive theory are scale invariant precisely at the 
ultraviolet nonattractive fixed point Unoo- The same remark was already pointed out in the 
m-axial Lifshitz universality SlJ following closely the behavior of the pure 0^ noncompeting 



case argument of BLZ 



301 ]. The several independent mass scales implementing independent 



renormalization group transformations in each type of competing subspace correspond to the 
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several correlation lengths appearing in each modulated ordered phase in the CECI model. 

We can consider direct applications of the present formalism in the first place. We can use 
the method to compute critical amplitude ratios of various thermodynamic potentials (as 
well as other universal observables constructed from the renormalized IPI vertex functions) 
above an below the Lifshitz temperature. Some uniaxial m-fold Lifshitz critical amplitudes 



have already been calculated at one- loop level [10|, |47|]. Even in this simpler case, not all 
amplitudes were computed and much more can be done. The motivation is to compare the 
results with new materials which might present Lifshitz points of generic higher character 
involving several independent length scales. Indeed, the comparison of the specific heat 
amplitude ratio of Ref.0] with the one obtained experimentally for the material MnP 



48l | yielded a remarkable agreement between the two results. Our hope is that the present 



method can help to unveil all the critical amplitude ratios, therefore generalizing the analysis 
for the m-axial Lifshitz with just two independent length scales. Moreover, we hope that 
numerical calculations such as those studied in high temperature series for the uniaxial 
Lifshitz case at fixed values of the space dimension and number of components of the order 



parameter [49[ shall be put forth in order to improve our understanding of the universality 
classes of the most general competing systems. 

Beyond the problem of critical phenomena occurring in generic competing critical sys- 
tems, the mathematical apparatus developed in the present paper might be useful to attack 
the perturbative aspects of many quantum field theories with higher derivatives that have 
been studied recently. It has been pointed out that field theories with higher space deriva- 
tives and some sort of mechanism like the Lifshitz condition which suppresses the second 
space derivatives eliminate the lack of unitarity and produces a well behaved theory without 



ghosts 



19|. It shares a similar renormalization group treatment: the time scale behaves in 



the same way as the noncompeting direction since it is quadratic in the derivatives, whereas 
the space directions have different scaling dimensions analogous to the Lifshitz competing 
directions. These Lorentz violating models have also been constructed in the cases of gauge 
;heories SOj with implications in the standard model and the analysis of neutrino masses 



5l|. 



Although not studied yet within our present investigation, the issue of quantum critical 
behaviors of Lifshitz points might also be worthwhile. An instigating aspect which appeared 
recently in the literature is the application of these ideas in examinating the quantum crit- 
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icality of membranes, which is analogous to the anisotropic m-axial n = 2 case with two 



length scales 



52| . This description corresponds to a new class of gravity models. It has been 



further speculated that the case n = 3 of third character anisotropic Lifshitz points can be 
studied in this framework and explains the basic features of quantum gravity at these Lif- 
shitz points [53]. This study adds more ingredients to understand the infrared modifications 
of gravity described by the ghost condensate previously introduced in Ref.[5J], which has 
a more transparent analogy with the Lifshitz CECI model: gravity could have attractive 
character at "short" distances, but might develop repulsive properties in the long distance 
limit. 

Furthermore, the conjectures put forward in 3l| can be further extended to the most 
general cases of arbitrary independent mass scales: each scale is characterized by an in- 
dependent Compton wavelength. This is kind of bizarre, since now more than one mass 
scale could characterize the corresponding "elementary particle" . This anisotropy could re- 
flect the way it is distributed over space(time). In other words, if space(time) would be 
anisotropic (as suggested by these many independent mass scales) this feature could leave 
a mark on particles propagating on it. The connection with the model discussed here can 
be easily viewed in considering a scalar quantum field with a Lifshitz-like condition which 
only keeps higher space derivatives and maintains solely second order time derivatives. A 
Wick rotation brings this theory to the statistical mechanics form of the Lagrangian density 
Eq.(2). In this "generalized Lifshitz space", the Laplacian operator is originally defined in 
d dimensions, but due to distinct types of competition axes (caused by forces with alternate 
signals) as realized in the CECI model, only (c? — X]n=2^") '^^ components (e.g., the time 
components) remain. The residual information that the field lived in (i-dimensions in the 
first place before dynamical effects break these pattern is encoded in the higher derivative 
terms in the Lagrangian density. 

Minkowski space corresponds to the limit of zero gravity, where quantum fields propagate 
in. Competition between repulsive and attractive components of gravity in the large distance 
limit could then produce a tiny effect, which could result in this Lifshitz space. When Wick 
rotated back, this would produce precisely another type of flat space limit where gravity is 
small but has an observable effect on it: Lorentz invar iance is broken in the bare propagator 
of the quantum fields, since now the timelike components of the momentum is quadratic as 
before, but the spacelike components have higher powers. Choosing independent metrics in 
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each competing subspace in the Lagrangian (2), the higher order term corresponding to the 
rrin dimensional competing subspace could be written as 5'rs(c^r)'^(<9s)", with r, s = 1, ...,n for 
Qrs = Srs- With this choice, the introduction of masses in dimensional reduction using the 
Siegel method would follow simply as discussed for each uniaxial subspace m„ = 1 without 
any further complication. We believe that the perturbative analysis of the examples above 
mentioned incorporating these new ideas might be worthwhile after the study described in 
the present work. 

Most of the remarks done in the anisotropic cases can be extended to isotropic points. 
From a realistic point of view, however, it is much more difficult to visualize examples 
in quantum field theories: as pointed out before, more than two time derivatives in the 
Lagrangian density leads to trouble with unitarity. We can make remarks concerning the 
method itself in comparison with the zero mass case previously discussed in Ref. 36|. 



The massive integrals computed using the orthogonal approximation are very simple 
and can be computed explicitly for generic number of neighbors coupled via competing 
interactions n. On the other hand, the exact massive integrals are rather involved, even at 
the one-loop order, within the proposal presented. Just as happens in the massless case, 
the one-loop massive diagram explicit computation with general n is identical to the result 
using the orthogonal approximation. The latter is exact at one loop. Deviations start at 
two-loop order and beyond. What is really remarkable using either the massless or the 
massive approach is that usual systems without competition are particular cases of generic 
isotropic criticalities with n = 1, i.e., they are Lifshitz behaviors of first character. 

In summary, we evaluated the critical indices of generic higher character Lifshitz points 
using massive fields along with normalization conditions defining the renormalized theories 
at zero external momenta with many independent masses responsible for independent renor- 
malization group transformations in each competing subspace. The results turn out to be 
identical to those obtained previously using a massless framework. Thus, universality is 
corroborated once again. 

XI. ACKNOWLEDGMENTS 

We would like to thank partial financial support by CNPq from Brazil. 



51 



APPENDIX A: FEYNMAN GRAPHS FOR ANISOTROPIC BEHAVIORS 



The required diagrams corresponding to one-, two- and three-loop integrals in momenta 
space shall be computed using dimensional regularization and Feynman parameters as our 
main tools in their solution. Solely the one-loop integral associated to the four-point IPI 
vertex part can be solved exactly for anisotropic behaviors. Multiloop graphs can be eval- 
uated through the generalized orthogonal approximation. The main step is actually quite 
simple: the loop momenta characterizing a certain competition subspace in a given bubble 
(subdiagram) do not mix to all loop momenta not belonging to that bubble. In other words, 
the loop momenta in a given subdiagram is orthogonal to all loop/external momenta ap- 
pearing in other subdiagrams. The normalization conditions defined in the text indicates 
that the integrals to be worked out should be determined at zero external momenta. The 
minimal set of Feynman integrals to be solved are presented with increasing order in the 
number of propagators, namely the one-loop integral contributing to the four-point function 



/ 



Et2(fc?n))"+(?)^ + l] 



(Al) 



the two- loop contribution to the two-point function I'^ — ^^^^^^|p=X(„)=o 
(= ^^^^^^|p=K(„)=o), where 73(P,X(„)) is the integral 



X r > (A2) 

the two-loop contribution to the coupling constant 

dd-i:^n=2 mnq^d'^-E^=^ "*"?2n^=2C^"^"^l{n)n^=2C^'""^2(n) 



-I 



(5f + Et2(%))" + l)" 



, (A3) 
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and finally the three-loop integral 1'^ = ^^^^^|p=X(„)=o(= ^^^^fi^\p=K(„.=o) with 

(n) 

I^{P,K(^n)) representing the graph 



' {qI + Et2(%))" + 1) (?i + Et2(%))" + 1) {4 + Et2(%))" + 1) 

[(gi + ?2 - P)^ + Et2((^l{n) + hin) - K^n)rT + 1] 

[(?! + 53 - P)^ + Et2((^l{n) + hin) - K^n)YT + 1] ' ^^^^^ 

We stress that in the expressions for and I^, P is the external momenta perpendicular 
to the several types of competing axes (,i-e., it belongs to the nii subspace) whereas K(^n} 
{n — 2,...,L) is the external momenta representing the nth m„-dimensional competing 
subspace. They are only needed to compute the derivative of these integrals, but shall be 
set to zero after that operation. 

One should keep in mind that whenever a loop integral is performed, the geometric angu- 

T c ran \ 

lar factor ['5(rf_j^^_^„„)r(2 - ^^^^2 i^)(n^=2 2« )] appears but can be omitted in the 
final result of the integrals by a redefinition of the coupling constant. Performing this redef- 
inition and applying the mathematical tools explained above, we learn that those divergent 
integrals have the following representation in terms of dimensional poles corresponding to 
their e^-expansion values: 

h = -[1 + {hm^ - l)eL] + 0(6^), (A5) 

7; = ^[1 + {2h^, - ^)6i] - 1/ + 0(6^), (A6) 
/4 = ^ (1 + (2 /i^, - l)e^ + 0{el)) , (A7) 
4 = ^[1 + i^hm, - \)eL + 0(ei] - (A8) 



where = 1 + ^» and 

ln[x{l — x)] 



r\ . 1 ln[x{l-: 
/ = / dx( z + 7 — 



xW' 



(A9) 



is the same integral occurring in the original work by BLZ and in the m-axial Lifshitz 
anisotropic integrals. This integral drops out in the calculation of the critical exponents by 
extensive cancellations in the renormalization functions at the fixed point. As pointed out in 
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the previous massive approach for the m-fold anisotropic Lifshitz criticahty, its appearing is 
a pecuharity of the orthogonal approximation. The normahzation conditions and our choice 
of the subtraction point are responsible for this nonuniversal function, since it does not show 
up in the e^-expansion of these integrals in the massless case. 



APPENDIX B: ISOTROPIC GRAPHS IN THE GENERALIZED ORTHOGO- 
NAL APPROXIMATION 

We shall pursue the analogy with the anisotropic case in order to compute the isotropic 
integrals {d — m„) in the orthogonal approximation. There is only one type of sub- 
space (/competition axes) to be integrated over and the parameter is now e„ = 4n — d. 
The minimal set of loop integrals now read: 

h = / o , Bl 

_ /■ d^"k,d^"k2 

' J (((A;i + A;2 + ir)2)- + l)p2)« + l]pi)« + l] ' ^ ^ 



d^^kid^'-ki 



[{kjY + 1]' {{klY + 1) {{{ki + k2fY + 1) ' 



(B3) 




+ A;2 + KfY + l) {{{k^ + k^ + KfY + l){{klY + l){{klY + l){{klY + 1) 



(B4) 

One should keep in mind that the derivatives with respect to of the two-point functions, 
namely l'^ and /g, are actually the objects of interest in our discussion. Recall that these 
integrals are evaluated at X = 0. Therefore, using the same technology as before, the 
calculation is even simpler than in the anisotropic case. Again, every time a loop integral is 
performed an angular factor takes place. The geometric angular factor to be absorbed in a 
redefinition of the coupling constant in the isotropic integrals is Smn ? which corresponds to 
the area of the unit sphere m. d — rUn dimensions. Omitting this factor, we encounter the 
following results for the required integrals: 

^2 = -[l-^en]+0(en), (B5) 
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where / is the same integral taking place in Appendix A. Similarly to the situation found 
in the anisotropic cases, it does not contribute to the isotropic critical exponents as well. 
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